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CONCEPTS FOR GENERATING OPTIMUM VERTICAL FLIGHT PROFILES 


John A. Sorensen 

Analytical Mechanics Associates, Inc. 
Mountain View, California 94043 


SUMMARY 


The objective of this project is to develop and evaluate one or more 
algorithms and flight management concepts for the on-board minimization 
of fuel and/or direct operating costs. These concepts are to be used for 
steering a CTOL aircraft in the vertical plane between fixed origin and 
destination airports along a given horizontal path. 

In this report, a summary is first made of work that has been con- 
ducted to date concerning the generation of optimum vertical profiles. 

Then, algorithms for generating optimum vertical profiles are derived and 
examined. These algorithms form the basis for the design of on-board 
flight management concepts. The variations in the optimum vertical pro- 
files (resulting from these concepts) due to variations in wind, 
takeoff weight, and range-to-destination are next presented. Further con- 
siderations for mechanizing two different on-board methods of computing 
near-optimum flight profiles are then outlined. Finally, the results are 
summarized, and recommendations are made for further work. 

Two appendices are Included which give technical details of optimum 
trajectory design, steering requirements for following these trajectories, 
and off-line computer programs for testing the concepts. The two computer 
programs developed are called OPTIM and TRAGEN, and they are available from 
the Computer Software Management Information Center, Barrow Hall, University 
of Georgia, Athens, Georgia 30601. 


V 




TABLE OF CONTENTS 


Page 


I INTRODUCTION 1 

II BACKGROUND 5 

III COMPUTATION OF THE NOMINAL OPTIMUM VERTICAL PROFILE ... 15 

Derivation of the Optimization Algorithm 15 

Generation of Nominal Optimum Trajectories 19 

Example Optimum Trajectories 23 

IV PROGRAM VERIFICATION AND SENSITIVITY RESULTS 31 

Optimization Verification 32 

Sensitivity Analysis . 37 

V MECHANIZATION OPTIONS 49 

Off-line Optimum Profile Generation 49 

On-board Optimum Profile Generation 51 

Method Comparison 53 

VI SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS 55 

Summary and Conclusions 55 

Recommendations 57 

APPENDIX A. TRAJECTORY OPTIMIZATION USING THE ENERGY STATE 

METHOD 63 

APPENDIX B. AIRCRAFT EQUATIONS OF MOTION AND AUTOPILOT MODELS 87 

REFERENCES 97 


yii 




LIST OF FIGURES 


Figure Page 

1. Projected Free World Petroleum Production Capacity vs Consumption. 2 

2. Free World Air Carrier Jet Fuel Consumption Projection 2 

3. Contours of Constant Energy Increase per lb. of Fuel Burned 

and Minimum-Fuel Energy-Climb Path 8 

4. One Solution to the Minimum Fuel-Fixed Range Problem 10 

5. Time-Fuel Tradeoff for Fixed Ranges 11 

6. Minimum Cruise Fuel Consumption per Unit Range for the DC-10 ... 12 

7. State Variables for an Optimum Climb Profile of a 200 n.mi. 

Range Flight 24 . 

8. State Variables for an Optimum Descent Profile of a 200 n.mi. 

Range Flight 26 

9. Comparison of Minimum Direct-Operating-Cost Flight Profiles. ... 29 

10. Effect of Range Variation on Optimum Climb Flight Path Angle 

and Mach Number 39 

11. Effect of Initial Weight Variation on Optimum Climb Flight 

Path Angle and Mach Number 40 

12. Wind Profiles Used to Study Sensitivity Effects on Optimum 

Profiles 42 

13. Effect of Wind Variations on Optimum Climb Flight Path Angle 

and Mach number 43 

14. Effect of Wind Variations on Optimum Descent Altitude and 

Mach Number as Functions of Range-to-go 44 

15. Block Diagram Depicting On-Board Climb Steering Based on 

Computing the Nominal Optimum Profile Off-line 50 

16. Block Diagram Depicting On-Board Climb Steering Based on 

Computing the Nominal Optimum Profile in the 

Background 52 

A.l Assumed Structure of Optimum Trajectories 67 

A. 2 Measurement of Range from the Origin or to the Destination .... 69 

A. 3 Optimum Cruise Cost as a Function of Cruise Energy 73 

A. 4 Optimum Profile Energy vs Range for > 0 at Cruise 75 

A. 5 Optimum Profile Energy vs Range for = 0 at Cruise 75 

A. 6 Plot of Thrust and Drag vs Airspeed at a Particular Altitude ... 79 

A. 7 Optimum Cruise Cost as a Function of Altitude for Cruise 

Weight of 54,431 kg (120,000 lb) 79 

A, 8 Optimum Cruise Cost and Cruis.e Altitude as Functions of 

Cruise Weight for the 727 Flying into a Particular Head 

Wind 

ix 



List of Figures, Continued 


Figure Page 

A. 9 Relationship Between the Cruise Cost Parameter ij; and the 

Associated Range of Flight 86 

B. l Vector Diagram of Longitudinal Forces 88 

B.2 Elements of the Longitudinal Aircraft Model 91 

B.3 Control Loops for Perturbation Control of Airspeed and 

Flight Path Angle 93 


X 



LIST OF TABLES 


Table Page 

1. Comparison of Flight Costs for 200 n.mi. Range with 

Varying Wind Conditions 30 

2. Comparison of Optimization and Trajectory Generation 

Program Climb Results for Different Range Flights 35 

3. Comparison of Optimization and Trajectory Generation 

Program Climb Results in the Presence of Winds 35 

4. Comparison of Optimization and Trajectory Generation 

Program Results Descending from 10668 m (35000 ft) 

to 3048 m (IQOOO ft) 36 

5. Comparison of Optimization and Trajectory Generation 

Program Results Descending from 10668 m (35000 ft) 

to 1.8 m (6 ft) 37 

6. Initial Weight Error Cases Considered - Climb Profile .... 46 

7. Wind Profile Error Cases Considered - Climb Profile 47 




SYMBOLS 


Roman Symbols 


C,C. 


D 


t 

D,D 

c 

DOC 

d 


d ,d 
up n 

E 

E ^ 
copt 


^up’^c’^dn 


g 


unit direct operating cost 
drag coefficient 
cost of fuel C$/lb) 
lift coefficient 
cost of time ($/hr) 
drag 

direct operating cost combining fuel and time costs 
cruise distance 

total distance traveled 

climb and descent range 

specific energy 

optimum value of specific energy during cruise at the 
optimum altitude 

fuel burned 

fuel burned during climb, cruise and descent 

constants in the transfer functions relating perturbations 
in Y to perturbations in a 

acceleration due to gravity 

the Hamiltonian used for optimization 


H ,H, 
up dn 

h 

J 

K^,K2,K3,K^- 

L 


value of Hamiltonian during climb and descent 
altitude 

cost function to be minimized 
control law gains 

lift 


L(x,u,t) 

m 

P 

Pl’^2 


cost per unit time along a trajectory 
aircraft mass 

ratio of cruise cost ijj to optimum cruise cost ip^ 
origin and destination values of range 


opt 


xlii 



q 

- 

pitch rate 

R . ,R 
min max 

- 

values of range where P = 1.5 and 1.01 on optimum trajectories 

r 

- 

ratio of fuel to time costs (C-/C. )/(l+C„/C. ) 

r t K t 


- 

range-to-go 

s 

- 

Laplace transform operator 

T 

- 

thrust 

t 

- 

time 


- 

time at beginning and ending of cruise 


- 

initial and final times of trajectory 

v,v ,v 

a’ c 

- 

airspeed 

V .V, 

up dn 

- 

airspeed during climb and descent 

V 

w 

- 

longitudinal component of wind 

V .V , 
wup wdn 

- 

wind speed during climb and descent 

W,W. ,W^ 
i f 

- 

weight of the aircraft; initial and final weight 

W 

c 

- 

weight during cruise 

• 

w 

- 

fuel flow rate 

X 

- 

range traveled during climb and descent 

X ,x , 

up dn 

— 

components of range traveled during climb and descent. 



Greek Sjnnbols 

a 

- 

angle-of-attack 

Y.Yc 

- 

flight path angle (with respect to air mass) 


- 

inertial flight path angle 


- 

throttle setting 

A 

- 

adjoint variable or costate used for optimization 

7T 

- 

throttle position or thrust control 

TT ,7T , 

up dn 

- 

throttle setting during climb and descent 

C 

- 

damping term in transfer function 

T 

- 

tlme-to-go 


- 

constants in the transfer functions relating perturbations 
in Y and to perturbations in a 


xiv 



terminal cost function 

cost per unit distance of cruise 

natural frequency in transfer function 
specific fuel consumption (Ib/lb thrust) 


XV 




I 


INTRODUCTION 


Virtually every segment of civilization is being affected by the 
energy crisis, and we are awakening to the fact that we live in a system 
with finite resources. Our natural supplies of hydrocarbon fuels are 
being rapidly consumed, so that it has become mandatory that we seek to 
develop every feasible means to conserve fuel. The energy crisis and 
its effect on developing technology is especially impacting the air trans- 
portation industry which accounts for about 4 percent of our national 
petroleum consumption [1]. 

The current motivating factor for our air transportation industry 
to conserve fuel is the escalating cost of jet fuel. To the commercial 
air carrier, reduction in fuel usage means reduction in operating costs 
and an increased profit margin [2]. For military air operations, fuel 
reduction also implies increased aircraft range, greater payloads, more 
missions, and an enhanced defense reserve [3,4]. The U.S. Air Force has 
an annual fuel usage of over 4 billion gallons [3], while our domestic 
commercial fleet used about 10 billion gallons in 1978 [1,2]. 

Figures 1 and 2, taken from Ref. 2, illustrate the dilemma that 
exists for commercial aviation - increasing air transportation demand 
versus decreasing production capacity. Figure 2 also illustrates that 
accelerated technological advancements can potentially reduce future 
usage of fuel for air transportation by more than 30%. One of these 
advancements is developing the ability to fly each aircraft along pro- 
files which minimize total fuel consumption. 

Commercial airline flight profiles are presently generated using a 
pseudo-optimization procedure. Given a reasonable knowledge of the jet 
stream, computer programs based on dynamic programming techniques test 
a series of horizontal paths between the departure and destination points 
to find the one providing minimum time or fuel usage. The horizontal 
paths tested are predetermined, and they typically consist of a series of 
Interconnected segments between VOR/DME and TACAN stations. The vertical 
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Figure 2. Free i7orld Air Carrier Jet Fuel Consumption Projection [2]. 

2 




path then consists of a fixed Mach number/indicated airspeed profile for 
climb and descent combined with a series of step changes in altitude 
during cruise. The altitude changes are governed by predetermined fuel 
burnoff rates and ATC considerations. This report is focused on various 
aspects of computing Improved vertical profiles. 

There is a pressing need to reduce fuel usage and costs for all 
types of flights. In addition, increasing labor costs and penalties on 
air travel time make it desirable to minimize the flight time. These 
factors motivate the search for a flexible procedure for generating a 
vertical reference profile which can minimize the total cost of an 
operation between two airports. 

This procedure should be capable of minimizing fuel usage, flight 
time, or a combination of these variables (direct operating cost - DOC). 
Furthermore, it should be computationally efficient and implemented so 
that it can be used by the pilot while airborne to change the flight pro- 
file in case of change in weather conditions or final destination. 
Finally, this procedure should serve as a means of automatically driving 
the aircraft control surfaces and engine throttle settings. The work 
summarized in this report has been directed at developing such a pro- 
cedure that can be mechanized on-board to achieve the optimum path. 

This report is organized as follows; 

a) A summary is made of publicly available work that has been 
conducted to date concerning the generation of optimum 
vertical profiles. 

b) Algorithms of methods for generating optimum vertical profiles 
are derived. These procedures form the basis for specifying 
on-board mechanization requirements. 

c) Sensitivity results are presented of the variation in the 
optimum vertical profile as a function of variations in wind, 
takeoff weight, and range-to-destlnation. These sensitivity 
results also impact the mechanization requirements. 
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d) Considerations are outlined concerning the on-board computation 
of the near-optimum profiles and the interface with the auto- 
pllot/autothrottle. This is the first step in the design of 

an experimental prototype system. 

e) Recommendations are made regarding further work that is re- 
quired so that the algorithms for minimizing DOC are available 
to all potential aircraft users and avionics designers. 

These items are covered in turn in the next five chapters of this report. 
Also, Appendices A and B are presented to give technical details of 
trajectory optimization, steering requirements for Implementing these 
trajectories,, and off-line computer programs for testing the concepts. 
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II 


BACKGROUND 


The basic principles of trajectory optimization are an outgrowth of 
applying the calculus of variations. These principles are outlined in 
the first part of Appendix A, and they are referred to throughout this 
chapter. The purpose of this chapter is to outline the known work that 
has been accomplished to this time using trajectory optimization methods 
for minimizing various cost functions including aircraft direct operating 
costs. This background serves to put into perspective the methods for 
generating near optimum trajectories described in Chapters III and IV. 

As mentioned in the Introduction, the first problem solved by the 
flight planner 1s to determine what horizontal path to fly between origin 
and destination to take best advantage of the winds aloft (jet stream). 
Here, they may be wishing to minimize fuel, time, or a combination of 
fuel and time (i.e., DOC). This problem is usually solved by assuming that 
a single flight level will be used and a particular wind vector field 
exists over this horizontal plane. This type of problem is typically 
solved using a computational technique called dynamic programming [7]. 
Computerized flight planning based on dynamic programming is available 
on a commercial basis to airlines and other users [8], Some large airlines 
currently have their own facilities for weather forecasting, prediction of 
wind fields, and planning of horizontal route structures. 

However, the primary subject addressed here is how to fly the aircraft 
in the vertical plane once the horizontal path has been specified. That 
is, the origin, destination, horizontal route, starting and ending 
altitudes and airspeeds, and takeoff weight are fixed. The problem then 
is to choose the vertical flight profile (altitude, airspeed, thrust, and 
flight path angle time histories) to minimize cost. 
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A set of non-linear differential equations describing a point mass 
model of the aircraft longitudinal motion is adequate to describe the ver- 
tical motion of the aircraft for optimization purposes [9]. This model 
is, for no wind, 


mV 

= 

T 

cos 

a 

- D 

- mg 

sin Y 

mVy 

= 

T 

sin 

a 

+ L 

- mg 

cos Y 

• 

h 

= 

V 

sin 

Y, 




• 

X 

= 

V 

cos 

Y, 




• 



• 





m 

= 

- 

w/g 

. 





( 1 ) 


Here, the state variables are the airspeed V, flight path angle y, altitude 
h, horizontal range x, and mass m. The control variables are the thrust T 
and angle-of-attack a which are both amplitude limited. The variables lift 
L and drag D are functions of h,V, and a. The fuel flow rate w is a function 
of h, V, and T. These equations are analogous to Eqs. (A.l) in Appendix A. 
The cost function to be minimized is of the form, 


J = f *(Cj»+C^)dt . 
t 

O 


( 2 ) 


where and are the unit costs of fuel ($/lb) and time ($/hr) . The 
final time t^ is free. The problem is to choose the sequence of controls 
(T and a ) that drive Eqs. (1), satisfy the initial and final constraints, 
and minimize the cost of the flight, as governed by Eq. (2). 

In addition to dynamic programming, other iterative procedures such 
as steepest descent, first and second order gradient, quasi-linearization, 
and neighboring extremal methods [5] have been developed and applied to 
solve aircraft trajectory optimization problems. All of these are com- 
putationally time consuming [3,9,10] and require a certain amount of art 
(gained by experience) in solving for the sequence of optimal controls. 
These problems led researchers to seek simplifications and assumptions to 
allow solving for an approximate optimum solution more quickly. 
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An early attempt at obtaining simplified solutions to the minimum time- 
to-climb and minimum fuel-to-climb problems was made by Rutowski [11] using 
graphical methods. He also suggested that the aircraft could adequately be 
represented in such a problem by its specific energy state, 

E = h + V^/2g . (3) 

The use. of this state assumes that kinetic and potential energy can be 
instantaneously interchanged. 


Bryson, Desai, and Hoffman [12] used the energy state approximation 
and differentiated Eq. (3) to obtain 


E = V(T - D)/mg 


(4) 


as the single equation of motion. This assumed that lift L balanced 
weight mg, and airspeed V was treated as a control variable. With these 
assumptions, simple solutions to the minimum time-to-cllmb problem were 
found by maximizing energy rate E, or 


max {V(T(E,V) - D(E,V))} 

V 

for all values of E. Minimum fuel-to-cllmb solutions were found by 
maximizing 


max {V(T(E,V) - D(E,V))/(-m)} 
V 


at each energy level. A typical solution to Eq. (6) is shown in Fig. 3. 
Solutions to maximum range for a given amount of fuel and maximum range 
for a given amount of time were also found. 

The work was extended to treat the minimum-fuel fixed-range problem 
by Schultz and Zagalsky [13-15]. In Ref, 13, they used range as the in- 
dependent variable and the state equations 
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MACH NUMBER 

Figure 3. Contours of Constant Energy Increase per lb, of Fuel 
Burned and Minimum-Fuel Energy-Climb Path [12] . 


dE/dx = (T(E,V,tt) - D(E,V))/mg, (7) 

dm/dx = -a(E,V,Tr)T(E,V,7r)/V, 

where 7T is the throttle position and o is the specific fuel consumption. With 
77 and V as control variables, they formed the Hamiltonian (See Appendix A), 

H = -aT/V + A(T-D)/mg. (8) 

Here, A is the adjoint (or costate) variable. This Hamiltonian is independent 
of X (and therefore constant) which eliminates the need to solve for A. H 
is maximized with respect to V and 77 and is used to show that for sufficient 
range, suboptimal trajectories consist of minimum fuel energy climbs, clas- 
sical steady state minimum fuel cruises, and maximum range descending glides. 
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In Ref. 14, Schultz and Zagalsky critiqued the energy state approxima- 
tion solutions by posing second, third, fourth, fifth, and seventh order 
models of the longitudinal aircraft equations of motion and then found and 
compared solutions to the associated fixed endpoint flight-path optimization 
problems. These included (a) minimum fuel-fixed range, (b) minimum time- 
fixed range, and, (c) minimum fuel-fixed range-fixed time trajectories. 

For the climb-cruise-descent trajectory, the steady cruise conditions were 
found to minimize. 


min (a D/V) = a 


D /V 

c c 


(9) 


Then, the climb solution was found by maximizing 


max (T - P)V/m 

h aT - a D V/V 

c c c. 

for TT and constant values of E. The descent solution minimized 
max 

min (D/m) 
h 


( 10 ) 


( 11 ) 


to maximize range by assuming thrust and fuel flow are zero. A typical 
solution is shown in Fig. 4. In Ref. 15, range was again used as the in- 
dependent variable, and the function 


]L _ max 

A ir,V 

derived from Eq. (8), was suggested as a means to formulate on-board 
algorithms for minimum time and fuel trajectories. 

Of corollary interest to the above work is the technical argument 
between Schultz [16] and Speyer [17] concerning whether steady-state 
cruise can ever be part of the optimum trajectory solution. However, 
from a practical point-of-view, cruise should be defined in a simple way 
(such as the steady state condition) so that on-board mechanization of 
the reference trajectory is facilitiated. 


(T - D)/mg 
H + w 


( 12 ) 
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MACH NUMBER 

Figure 4. One Solution to the Minimum Fuel-Fixed Range Problem [14]. 

Also of related interest are the efforts of Kelley [18], Calise [4,19], 
and Ardema [20] in applying singular perturbation theory to find rapid 
solutions to the "exact" trajectory optimization problems and to Improve 
the results of energy state methods. In Ref. [20], Ardema compares energy 
state, zeroth order and first order singular perturbation, and steepest 
descent solutions to the minimum time-to-climb problem. The first-order 
singular perturbation solution produced essentially the same trajectory 
but required about l/40th the computation cost of the steepest descent 
solution. Aggarwal and Calise [4] found zeroth order reference solutions 
for the minimum cost climb-cruise-descent problem for five different Air 
Force bombers and tankers flying over fixed ranges. 

Recent contributions to the practical solution of minimum cost vertical 
profiles have been obtained by Erzberger and his collegiies at NASA Ames 
Research Center [21,22,24,27]. Erzberger uses the assumption that aircraft 
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energy is monotonically increasing during climb and monotonically de- 
creasing during descent. Thus, E is chosen to serve as the independent 
variable. This result produces a constant value for the adjoint variable 
X on an optimum trajectory - that being the steady-state cruise cost. The 
result is minimization of the Hamiltonian at increasing energy values 
during climb, or 


min 

1T,V 


C. w + C - 1(V + V ) 
f t w 

(T-D)V/mg 


(13) 


where the controls are airspeed and throttle setting. Here, is the 
longitudinal component of the wind. ' (More details are given about the 
derivation of Eq. (13) in Chapter III.) The descent portion of the trajec- 
tory is generated backwards in time with energy rate constrained to be 
negative. As can be seen, Eq. (13) is the inverse of Eq. (10) when 
allowance is made for the cost of time, the longitudinal effect of wind, 
and the difference in controls. The algorithm developed from Eq. (13) 
was used to determine the tradeoff between fuel and time for fixed range 
flights as shown in Fig. 5. Equation (13) also can be incorporated as 
an on-board algorithm, as discussed in Chapter V. 



Figure 5. Time-Fuel Tradeoff for Fixed Ranges [21]. 
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It was further recognized by Erzberger (24) that optimum trajectories 
with thrust constrained to the maximum value for climb and to the idle 
value for descent have different characteristic shapes than those where 
climb and descent thrust are free. These characteristic shapes were 
derived by analysis of the transversality conditions, and they are sketched 
in Figs. A. 4 and A. 5 of Appendix A. Using variable thrust can reduce fuel 
consumption by about 1% for certain engines (where fuel flow rate is not 
linearly proportional to thrust) . It was also revealed by Erzberger and 
others [22,23] that with a shear headwind, there could be more than one 
optimum cruise altitude for the same range. This is illustrated in Fig. 6. 
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Figure 6. 



Bochem and Mossman [25] have investigated the implementation of flight 
path and thrust management strategies outlined in Ref. 24. This effort 
involved real-time simulation of the flight path guidance algorithms on 
a DC-10 simulator, and evaluation of the system performance in a simulated 
environment . 

The Boeing Company [26] is also studying the application of the above 
methods to generate near-optimum vertical profiles. One Boeing algorithm 
used to find reference trajectories which reduces fuel usage finds 


max (T-D) /mg , 

V (aT/V) - a D /V ’ 
c c c 


for climb and 


IT = ir 


max 


(14) 


min (D-T)/mg . 

V (a D /V ) - OT/V ’ 
c c c 


TT = 7T 


min 


(15) 


for descent. These can be seen as variations of Eq. (10) and a special 
case of Eq. (13) . 


Current work being conducted by NASA Langley Research Center in- 
cludes developing the on-board capability for making idle thrust fuel 
efficient descents to a metering fix at a desired time previously 
specified by the ATC system [27]. They have also studied advanced air 
transportation concepts using flying wings which may reduce fuel usage 
up to 90% [28]. NASA Ames Research Center is currently flight testing 
minimum noise/fuel descents for powered lift STOL aircraft [29]. 


In summary, there has been considerable research work that has been 
directed at the problem of reducing in-flight direct operating costs. 

Several companies today advertise that they sell such systems (black boxes) 
which can be used to reduce flight fuel consumption. However, the following 
details are normally missing in discriptions of these systems: what specific 
algorithm is implemented, how is it interfaced to sensors and other avionics 
equipment, what computer capacity is required, and what amount of fuel 
(and time) is actually saved by the system compared to a true optimum. It 
is the intent of this work to begin to provide such details for a near- 
optimum system. 
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Ill 


COMPUTATION OF THE NOMINAL OPTIMUM VERTICAL PROFILE 


In this chapter, an algorithm is derived which serves as the basis 
for generating an optimum vertical profile. This profile can either be 
generated before flight or as a background computation while in flight. 
Considerations of the mechanization are discussed in Chapter V. 

Derivation of the Optimization Algorithm 
Begin with the cost function given as Eq. (2) and repeated here. 

r *"f 

J = / (C^ w + C^)dt. (16) 

0 

That is, over the length of flight, it is desired to minimize a combina- 
tion of both fuel and time. These quantities cost C^ ($/lb) and C^ ($/hr) 
This cost function is general in the sense that if C^ is zero, the 
trajectory will be in minimum time. If C^. is zero, the trajectory will 
be minimum fuel. 

For general trajectory shaping, the flight path angle dynamics of Eqs 
(1) can be neglected, and it can be assumed that the angle-of-attack a is 
small. Then, the equations of motion of the aircraft in the longitudinal 
plane are 

V = (T-D)/m - g sin Y , (17) 

h = V sin Y , 

X = V cos Y + V , 

w 

where the longitudinal component of the wind velocity is included . 
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Now, based on previous researchers' experience, the specific energy 
and its derivative are defined according to Eqs. (3) and (4). The first 
and second of Eqs. (17) are replaced with 

E = V(T-D)/mg . (18) 

Consider again the cost function of Eq. (16). What is actually 
desired is to minimize the cost in flying from origin to destination P2, 
where the time of flight t^ is considered to be free. The cost function 
can be changed to have range as the independent variable by dividing Eq. 
(16) by the expression for range rate, or 


J 



(Cj w + Cj.) 

7; • dx 

V cos y + V 
' w 




The remaining state equation (Eq. (18)) also is changed as 


(19) 


dE ^ ^ ^ 

dx dt dt * 


= (T-D) /_V 

mg \ V cos Y + 


( 20 ) 


This neglects the rate of fuel flow as a function of distance. Practical 
experience [24] has shown that as long as the current value of aircraft 
weight (mg) is used in Eq. (20), this simplification is appropriate. 


Equations (19) and (20) are combined, using the procedures outlined 
in Appendix A to form the Hamiltonian, 


H 


dc ^ , dE 


w + + XV(T-D)/mg 

i. c 

V cos Y + V 
. ' w 


( 21 ) 
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As can be seen, this is a more elaborate Hamiltonian than assumed by 
Schultz and Zagalsky in Eq. (8). It allows for penalizing of time, it 
assumes no particular form for the fuel flow rate w, and it takes into 
account the horizontal component of the wind V^. 

Note, again however, that none of the terms in Eq. (21) are explicit 
functions of the independent variable x. Then, according to the arguments 
associated with Eq. (A. 12) in Appendix A, this Hamiltonian is a constant. 

This greatly facilitates obtaining the optimum profile solution. 

Again, assume that the thrust T is a function of the engine pressure 
ratio (EPR) setting Tr. The two control variables can then be considered to 
be TT and airspeed V for this simplified d 3 mamic model. (Later, these 
arguments will be extended so that the conventional controls - throttle 
and elevator deflection - are used). Both it and V have upper and lower 
bounds. By following Eq. (A. 11), the Hamiltonian H is minimized according to 

> 0 , (22a) 

> 0 , (22b) 

for all admissible values of it and V. For an assumed constant value for 
H, Eqs. (21) and (22) constitute three equations with three unknowns - 
TT, V, and A at every point along the optimum trajectory. 

Using the arguments of Ref. (21) - (24), assume that the vertical 
flight profile consists of three phases: 

(a) a climb portion where both altitude and energy are 
monotonically increasing with range ; 

(b) an essentially equilibrium cruise portion where thrust 
equals drag and lift equals weight; and 

(c) a descent portion where both altitude and energy are 
monotonically decreasing as range-to-go decreases. 






17 



For the cruise portion, the Hamiltonian can then be written as, 


H 

c 


Cj w + C 
V + V 


t 


w 


(23) 


For a given cruise altitude and weight, the EPR setting and airspeed are 
found by numerical search techniques such that the value of this expres- 
sion is minimized. This value of the Hamiltonian is then the constant 
value used to climb to and descend from this altitude. It also represents 
the minimum cost per unit distance traveled along the cruise path. 

(This implies that the cruise altitude and weight are known. These issues 
are also discussed later.) 


The constant value of H^ can be substituted into Eq. (21) , and Eq. 
(21) can be solved for the negative value of the unknown costate: 


- X = 


C. w + C - H (Vcos Y + V ) 
f t c w 

V(T-D)/mg 


(24) 


Thus, to minimize the cost function expressed by Eq. (21) requires that 
values of 7r and V be chosen such that Eq. (24) is minimized at every point 
along the climb and descent path. Equation (24) is the key algorithm to 
both offline and on-board computation of the optimum trajectory that mini- 
mizes direct operating cost. It can be seen again that this is a more 
elaborate version of the algorithm suggested by Zagalsky in Eq. (12). 

The work of Erzberger [21,22,24,26,30] has been based on an alternate 
approach to this problem. He used specific energy E as the independent 
variable with the assumption it always increases monotonlcally during climb 
and that it decreases monotonlcally during descent. Thus, the remaining 
state variable is the last of Eqs. (17), and Eq. (21) is replaced by 


H = ^ + 


dc 

dE 


dx 
^ dE 


w + Cj. + X (V cosY + VJ 
V(T-D)/mg 


(25) 


This is a repeat of Eq. (13) with a sign change on X. 
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Equation (25) can be compared to Eq. (24), and it can be seen that 
the Hamiltonian H and the costate variable X have reversed roles . Thus , 
these two solutions are duals of one another. For energy as the state 
variable, Erzberger shows, using transversality conditions, that X is con- 
stant and is equal to the negative value of the cruise cost. In this 
case, the Hamiltonian varies in magnitude because variables within Eq. 

(25) vary with the magnitude of the specific energy (the independent 
variable) . 


Generation of Nominal Optimum Trajectories 

Because Eq. (24) and (25) are essentially identical in character and 
concept, it was convenient to make use of the computer program previously 
developed by Erzberger and Lee [30] to procede with generating optimum pro- 
files. This program is based on the use of Eq. (25). To avoid confusion, 
the following description refers to minimizing the Hamiltonian of Eq. (25) 
with the costate variable X treated as the constant. The computer program 
uses principles described in Ref. 30 and is summarized in the second part 
of Appendix A. This program (referred to as OPTIM) which generates an 
optimum profile between origin and destination points was extended to have 
additional desired features, and its specific characteristics (from a pro- 
gramming point of view) have been documented in a companion users* guide 
[31] to this report. It is instructive here, to present in general terms, 
the elements and steps required to generate an optimum trajectory. 

The following quantities must be input to the program for generating 
points on the optimum profile: 

a) The aircraft initial takeoff altitude, airspeed, and weight, 

b) The final desired aircraft altitude and airspeed, 

c) The range between and heading to be followed by the aircraft 
from the origin to the destination points, 

d) The wind velocity and temperature profiles (if other than 
normal) as functions of altitude, and 

e) The values of the constants and in the cost function. 

' ft 
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The program must have access to numerical data (either tables or poly- 
nomials) which produce values of lift and drag coefficients C, and C^, 

Li D 

thrust T, and fuel flow rate w as functions of altitude, airspeed, 
temperature, and EPR setting. Also, for the particular aircraft being 
used, the program must have the following two sets of empirical equations: 

a) . A method of guessing the approximate initial cruise weight 

based on initial takeoff weight and the values of fuel and 
time cost; and 

b) . A method of guessing the approximate landing weight based 

on the final cruise weight. 

Finally, the program must have an efficient procedure which can be 
used to find the minimum value of the Hamiltonian (Eq. (25)) at each 
energy level by varying one or more variables (here , it and V) over their 
admissible regions. For the OPTIM program, the Fibonacci search technique 
is used. 

It must be mentioned that for some cases, there is no cruise section 
of flight, as is explained in Appendix A. In this case, the constant 
value of the costate variable refers to what would be the cruise cost at 
the highest point on the trajectory where the profile transitions from 
climb to descent. 


The basic steps that the optimum profile generation program follows are 

1) . The minimum value of cruise cost (Eq. (23) and Eq. (36)) is 

evaluated for various cruise weights and with altitude varying 
from sea level to ceiling height. The results are tabulated 
in a "cruise table". A typical example of cruise cost as a 
function of altitude is shown in Fig. A. 7. 

2) . Based on initial weight and values of and C^, the amount 

of fuel required for the climb is approximated (See Eq. (A. 38)). 
This produces cruise weight which is used to determine the 
optimum cruise cost and altitude from the cruise table. The 
cruise cost (A in Eq. (25); - ’I^(E^) in Appendix A) is required 

for minimization of Eq. (25). Also, choosing the cruise cost 
from the cruise table also produces the associated maximum 
values of energy E^ where the climb portion of flight ends. 
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3). The next step is to minimize the Hamiltonian (Eq. (25)) at 
« each point along the climb trajectory. This begins at the 
computed initial energy and stops with the precomputed final 
climb energy. The numerical procedure precedes as follows: 

a) . A step is made in specific energy (E^ ~ ^n 1 ^ * 

where the step size is held constant. 

b) . The range of acceptable values of V and tt are com- 

puted for this energy value . 

c) . Over these ranges, the Hamiltonian (Eq. (25)) is minimized 

with Y assumed to be small (cos y = 1) . This, then 
specifies the values of airspeed, lift, drag, thrust, 
altitude, fuel flow rate, and energy rate (Eq. (18)). 

d) . The time required to make this step is computed as 


At = AE/E . (26) 

e) . The flight path angle is then computed as 

Y = sin"^ (Ah/At) /V. (27) 

f ) . The value of range is increased by Ax, where 

Ax = (V cos y + V (h))At. (28) 

w 

g) . The aircraft weight is decreased by Aw, where 

Aw = w At. (29) 


This procedure is repeated until the initial cruise energy 
level is reached. This produces an accurate measure of the 
fuel burned during climb and the range traveled during climb. 

4) . A final approach (or landing) weight is approximated by 

estimating the cruise range, cruise fuel, and other related 
quantities (See Eqs. (A. 42) - (A. 46)). This procedure also pro- 
duces an estimate of the final cruise altitude, energy, and 
cruise cost. 

5) . The procedures for generating the descent profile are now 

followed using the same sequence as for Step (3) but by going 
backwards in time and constraining the energy rate to be 
negative. This is continued until the estimated final cruise 
energy level is reached. This produces an accurate estimate of 
the fuel burned during descent and the range traveled during 
descent . 


21 



6) . The cruise distance is now recomputed by subtracting the climb 

and descent ranges from the total required. This cruise distance 
is used to compute the final cruise weight and the final approach 
(landing) weight, 

7) . With the revised final approach weight, the descent trajectory 

is recomputed using Step (5). This also provides refinement of 
the cruise distance, the total fuel burned, and time required 
to travel from origin to destination airports. 

The above description is somewhat simplified in that profiles of 
short range flights (e.g., less than 250 n.mi.) require special computa- 
tions. Figures A. 4 and A. 5 in Appendix A depict the shapes of the 
trajectories for short range for two types of control strategies, If 
both V and IT are used as controls, the trajectories have the shapes de- 
picted in Fig. A. 5. For this case, with short range there is no cruise 
segment, and an iterative process is required to obtain the constant value 
of A and the maximum energy state used to compute the climb and descent 
profiles. 


If only airspeed is used as the control variable, then thrust is 
set at maximum climb value for climb and idle value for descent. For the 
727 aircraft model used in this study, this increased the cost of the 
typical short-range profile about 1%. However, it simplifies mechaniza- 
tion and it reduces the time required to compute the profile. Trajectories 
for this case have the characteristics depicted in Fig. A. 4. That is, 
they have cruise segments, but these segments occur at altitude and energy 


levels lower than the optimal value E 

copt 


Again, the program must iterate 


on the correct value of A, and Eq (A. 48) in Appendix A is used to compute 


the corresponding length of the cruise segment. 


The OPTIM program has been constructed so that various options can 
be exercised: 

1) . No wind or arbitrary wind 

2) . Optimize with V and 77 or V only 

3) . Arbitrary cost terms and 
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A) . With or without a speed constraint of 250 kt below 30A8 ja 
(10,000 ft.) 

5). Start in cruise or after takeoff. For starting in cruise, 

the cruise altitude can be set to a fixed value, or the pro- 
gram will compute the best value associated with the given 
cruise weight. 

In addition to generating optimum nominal profiles for the 727 aircraft, 
OPTIM can be used to examine the sensitivities in trip cost due to varia- 
tions in wind, range, takeoff weight, and the cost terms and C^. 

Examples of these variations will be given in the next chapter. 

The computations that take place in OPTIM for off-line generation of 
optimum vertical profiles are also the same as would be used by an on- 
board flight management oystem. For the on-board mechanization, the flight 
engineer would enter initial and final altitude and airspeed, initial weight, 
range, and heading. Wind and temperature profiles and aircraft weight 
estimates would be updated during flight. The seven basic steps listed 
earlier would be recomputed every few minutes to provide a current reference 
trajectory to be followed. 


Example Optimum Trajectories 

Figures 7 and 8 show plots of various state variables as functions 
of time (or time-to-go) for optimum climb and descent of a 727 aircraft 
traveling 200 n.mi. in range. For this example, the indicated airspeed 
was constrained below 250 kt below 3048 m (10,000 ft). Also, Optimization 
was done by varying V only, with 7T set to the maximum value for climb and 
idle value for descent. The cost constants were $0. 265/kg ($0. 12/lb) 
fuel and $600. 00/hr. As can be seen, when 3048 m (10,000 ft) is reached, 
the aircraft levels off and accelerates before continuing to climb. The 
opposite is true for the descent profile. For these profiles, the reference 
flight path angle is rough because of the inherent assumption that either 
altitude, airspeed, or both can be abruptly changed to obtain a change in 
energy. Flight path angle can be smoothed after it is generated to pro- 
vide a more flyable path. Note also that no where on either the climb or 
descent is Mach number constant (contrary to profiles that are specified 
in aircraft handbooks and normally flown). 
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Figure 9 compares the altitude vs. range profiles for three cases of 
a 727 aircraft covering 200 n.mi. For these, both V and tt were varied. 
However, a head wind and tail wind were added in two of the cases to pro- 
duce the differences. The wind magnitude profile is also shown in Fig. 9, 
and it is representative of the Denver wind measured in August, 1977. 

Table 1 makes various range and cost comparisons of the three pro- 
files shown in Fig. 9, based on values of of $0. 138/kg ($0. 0628/lb) and 
of $500. 00/hr. It is seen that the head wind caused a total cost in- 
crease of $24.95 (or 4.3%) while the tail wind decreased the cost $23.59 
(or 4.1%) compared to the no wind optimum case. No procedures were avail- 
able for predicting the cost of these flights if the 727 aircraft were to 
follow the profile specified in the aircraft handbook. Data, such as 
presented in Table 1, are printed for each run of OPTIM so that the trajec- 
tory costs and other characteristics can be readily evaluated (Again, see 
Ref. 31). 
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Figure 9. Comparison of Minimum Direct-Operating-Cost Flight Profiles 


29 


Table 1. Comparison of Flight Costs for 200 n,mi. Range with 
Varying Wind Conditions 




Case 



No Wind 

Head Wind 

Tall Wind 

Fuel Used (kg (lb) ) 




Climb 

1945.8 

(4289.60) 

2013.0 

(4437.82) 

1805.2 

(3979.66) 

Descent 

272.2 
( 600.19) 

295.1 
( 650.57) 

284.9 
( 628.13) 

Total 

2218.0 

(4889.79) 

2308.1 

(5088.39) 

2090.1 

(4607.80) 

Range (n.mi.) 




Climb 

115.51 

115.90 

107.36 

Descent 

85.24 

85.78 

94.60 

Total 

200.75 

201.68 

201.96 

Time (min; sec) 




Climb 

17:31 

18:33 

15:41 

Descent 

15:05 

15:33 

16:11 

Total 

1 

32:36 

34:06 

31:53 

Cost ($) 




Climb 

414.51 

432.40 

379.91 

Descent 

163.29 

170.35 

174.30 

Total 

577.80 

602.75 

554.21 

Cost/Distance 
($/n.mi. ) 




Climb 

3.59 

3.73 

3.54 

Descent 

1.92 

1.99 

1.84 

Total 

2.88 

2.99 

2.74 
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PROGRAM VERIFICATION AND SENSITIVITY RESULTS 


The previous chapter and Appendix A describe an efficient way in which 
near’-optimum flight profiles can be generated without using time-consuming 
numerical techniques such as described in Refs. 9 and 10. Instead of 
iteratively solving the two-point boundary value problem, assumptions 
were made so that the dynamics are simplified to two state variables. One 
state variable (either range or energy) becomes the independent variable, 
and the other was included as the single state in the Hamiltonian. Thus, 
the problem of solving for minimum cost profiles reduced to algebraic 
minimization of the Hamiltonian at each point along the profile. 

Now, although this method is a convenient way of generating a trajec- 
tory, the resulting trajectory must be verified by using a more accurate 
model of the aircraft dynamics, Verification implies that: 

1) The reference trajectory that is generated must be flyable 
when the full aircraft equations of motion and constraints 
are taken into account, 

2) The trajectory cost (Eq. (2)) as predicted by OPTIM must be 
essentially Identical to that experienced by simulating more 
complete aircraft equations of motion . 

Thus, the first objective of this chapter is to describe a companion 
program to OPTIM which was developed for verification of the optimization 
results. This program is referred to as TRAGEN (for trajectory generation). 

With OPTIM and TRAGEN as computer tools, the user has the capability 
to study the characteristics of optimum profiles in great detail and to 
examine alternate ways these profiles can be implemented on-board. Trajec- 
tory characteristics are obtained by making sensitivity studies which is 
the subject of the second part of this chapter. 



Optimization Verification 


To achieve optimization verification required the development of a 
companion computer program to OPTIM which can be used to simulate the 
longitudinal trajectory of an aircraft commanded to follow the reference 
path output from OPTIM. Some details of this program (TRAGEN) are pre- 
sented in Appendix B, and a separate user^s guide for this program has 
also been written [32]. 


In addition to verification of the optimization program’s results, 
the TRAGEN program has the following utility; 

1) . It provides a means for testing guidance laws for steering 

the aircraft to follow the input reference trajectory. 

2) . It enables study of the effect of following an incorrect 

reference trajectory. For example, if the OPTIM results were 
based on one particular wind profile and initial aircraft weight, 
and a different weight and wind profile actually existed, the 
TRAGEN simulation would allow assessment of the effect of these 
errors on trajectory cost. 

3) . It can be used to determine the flight cost that would result 

from the aircraft being commanded to follow a reference trajec- 
tory suggested in the manufacturer's aircraft handbook. For 
example, for climb, handbook reference trajectories usually 
consist of following a constant indicated airspeed until a 
given Mach number is reached. Then, the reference trajectory 
follows this fixed Mach number until the reference cruise 
altitude is reached. 

4) . It can be used to test perturbation control schemes for removing 

the effect of wind gusts and other non-nominal performance 
sources (navigation errors, transient temperature profiles, 
non-standard engine performance) . 

5) . It is expandable to test candidate on-board mechanizations of 

a system of equations for generating the near-optimum vertical 
profile. 


A five state-variable model of the aircraft is currently used in 
TRAGEN to simulate longitudinal motion. State variables are altitude, 
altitude rate, longitudinal range, airspeed, and aircraft mass. Currently 
neglected are the rapid transient dynamics of throttle response, angle-of- 
attack, and pitch rate (6^, a, q). The throttle is assumed to be set so 
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that maximum thrust is achieved during climb and idle thrust is used 
during descent. The cruise phase is not simulated. The altitude control 
variable is taken to be the angle-of -attack which has maximum and minimum 
limits. This degree of sophistication is adequate for testing the OPTIM 
results . 

The TRAGEN program can readily be expanded to include throttle 
dynamics and short period dynamics. This would be required for further 
study of autothrottle and autopilot design to steer the aircraft to follow 
input reference trajectories. The control variables would be throttle 
position and elevator deflection for this expanded capability. A require- 
ment for implementing this expanded simulation would be to obtain the 
necessary stability and control derivatives to complete the dynamic model. 


The specification of reference profiles used in TRAGEN is based on 
using altitude as the independent variable for climb and range-to-go to the 
destination as the independent variable for descent. For climb, the 
reference trajectory consists of specifying airspeed and flight path angle 
(with respect to the air mass) as piecewise linear functions of altitude. 

At the 3048 m (10000 ft) point, the aircraft is commanded to level off and ac- 
celerate until the airspeed is reached where the climb should again continue. 

To generate the control law to follow the commanded climb profile, a 
linear perturbation model was made of the dynamic equations , 

T cos a - D(a,h,v) - W sin y = mV , (30) 

T sin a + L(a,h,v) - W cos y = mVy • 

The perturbation equations and transfer functions from Eqs. (30) are 
given in Eqs. (B.IO) and (B.ll) in Appendix B. Here, it is assumed that 
a perturbation 5a to the nominal angle-of-attack can be used to obtain 
the desired perturbations in flight path angle (6y) and airspeed (6V) to 
maintain the aircraft on the desired reference climb profile. 
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Because both commanded values of airspeed and flight path angle are 
assumed to vary linearly with altitude, they are ramp functions, and a 
Type 1 control system is suggested. The original climb control law was 
thus set to be 

K K 

6 a = (K^ + j?-) (V^ - V) + CK 3 + -^) (Y^ - Y) (31) 

where V and y are the commanded reference values of V and y. K, , K„, K_, 
c ' c '12 3 

3 ^>^d are control gains. The commanded values were derived according to 
Eqs. (B.12) and (B.13). 

Experience with this control law revealed two points : 

1) The integral control gain K 2 on airspeed error was not needed 

and did not particularly improve performance. Thus, it was 
nominally set to zero. 

2) Constant values of the gains K^, and could be selected 

to provide good, stable performance throughout the climb phase. 
Thus, there was no need to have altitude dependent gains pro- 
grammed , 

No attempt was made to select the control gains so that the perturba- 
tion response was optimized. The main objective was to obtain a set of 
gains which caused the aircraft simulation to track the input reference 
trajectory adequately (which was accomplished). Gain selection for control 
response optimization should remain as a task to be conducted when an actual 
autopilot/autothrottle is being implemented where the complete aircraft 
dynamics are considered. 

Several optimum trajectories were computed by using OPTIM and then 
subsequently used as inputs to drive the TRAGEN simulation. For example. 
Table 2 presents a comparison of OPTIM and TRAGEN results for the climb 
portion of a 727 initially weighing 61236 kg (135000 lb) , and traveling a 
total range of 150, 225, and 275 n.mi. with no wind. As can be seen the 
match is exceptional. The same degree of comparison was found using the 
two programs for different range and different initial weight flights. 
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Table 2. Comparison of Optimization and Trajectory Generation 
Program Climb Results for Different Range Flights. No Wind. 
61236 kg (135000 lb) Takeoff Weight. 


Range 

Fuel Burn - 

Time 

- Sec 

Altitude 

m 

‘ (ft) 

Airspeed 

m/s 

(£t/e) 

Range 

1 

s 

n. ml. 

OPriM 

TRAGEN 

OPTIM 


OPTih 



TRAGEtJ 

OPTIM 

TRAGEN 

150 

U05 

1405 

682 

680 

7597 

7597 

230 


125127 

125054 

(3097) 

(3097) 

(24924) 

(24926) 

(756) 

(752) 

(410522) 

(410281) 

225 

1931 

1931 

1060 

1060 

9760 

9760 

244 


211441 

212936 

(4258) 

(4258) 

(32021) 

(32021) 

(801) 

(800) 

(693704) 

(698609) 

275 

2075 

2073 

1185 

1178 

10319 

10321 

244 


242116 

241033 

(4574) 

(4570) 

1 

(33855) 

(33861) 

(802) 

(800) ! 
1 

(794345) 

(790790) 


The results with head and tail winds were not as close. Table 3 pre- 
sents a comparison of OPTIM and TRAGEN results when the wind profile of 
Fig. 9 was used as both a head and tail wind for the 225 n.mi. range flight. 

In both cases, the more detailed simulation from TRAGEN shows that it takes 
a longer time period (14-22 sec) and greater range (4000 - 5000 m (14000 - 
15000 ft)) than predicted by OPTIM (However, range traveled is not signifi- 
cant for climb.) The biggest discrepancy is the 2.6% extra fuel required 
for the tail wind case. In the future, the modeling simplifications of 
OPTIM and the steering accuracy of the TRAGEN control law should be investigated 
to resolve this point. 


Table 3. Comparison of Optimization and Trajectory Generation 
Program Climb Results in the Presence of Winds. 225 n.mi. 
range. 61236 kg (135000 lb) 


Type 


Time 

- Sec 




Range 

m 

~ (ft) 

Wind 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

Head 

2084 

(4594) 

2082 

(4591) 

1180 

1202 

10294 

(33772) 

10362 

(33996) 

244 

(802) 

244 

(802) 

225177 

(738770) 

229818 

(753997) 

Tall 

1701 

(3751) 

1747 

(3851) 

897 

911 

9145 

(30003) 

9074 

(29770) 

233 

(766) 

235 

(772) 

183238 

(601174) 

187433 

(614936) 
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For the descent portion of the flight, the more likely variable - range- 
to-go to the destination point - was used as the independent variable for 
specifying the reference path. In this case, the descent altitude was 
specified as a function of range. This defined the inertial flight path 
angle that the aircraft should be on which is defined by Eq. (B.16). In 
the current version of TRAGEN, no attempt was made to control airspeed 
during the descent. The control law to compute the perturbation to the 
angle-of-attack was of the form 

K, 

6a = (K3 + - Y^-) . (32) 

where Yj-^ is the commanded inertial flight path angle. 

Despite the lack of control of airspeed during descent, the descent 
steering provides the same degree of matching performance between OPTIM and 
TRAGEN results as was experienced during the climb simulation (where air- 
speed was controlled) down to 3048 m (10,000 ft). Table 4 compares OPTIM 
and TRAGEN results at about 3048 m (10,000 ft) for a 727 descending with 
idle thrust from 10668 m (35,000 ft). The three cases are with no wind and 
the Denver head and tail winds shown in Fig. 9. The biggest error is in 
airspeed where variations of +5.8 to -1.5 m/s (+19 to -5 ft/sec) are seen. 


Table 4. Comparison of Optimization and Trajectory Generation 
Program Results Descending from 10668 m (35000 ft) to 3048 m 
(10000 ft). 


Type 


Time-to- 

-go- 8 


BH 



Range-to 

m 

-8° - (ft) 

Wind 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

No Wind 

117 

117 

444 

447 

3154 

3148 

149 

151 

56704 

56598 

(259) 

(258) 

(10348) 

(10328) 

(490) 

(497) 

(186036) 

(185688) 

Head 

110 

110 

444 

449 

3154 

3148 

149 

148 

55465 

55494 

(242) 

(242) 

(10348) 

(10328) 

(490) 

(485) 

(181973) 

(182066) 

Tail 

123 

124 

440 

441 

3154 

3148 

149 

155 

57911 

57856 

(271) 

(273) 

(10348) 

(10329) 

(490) 

(509) 

(189998) 

(189816) 
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These same cases are compared again at the altitude of 1.8 m (6 ft) in 
Table 5. The head wind case has a -3.4 m/s (-11 ft/sec) difference in air- 
speed which causes 7 sec time and 4.5 kg (10 lb) fuel burn differences at 
this point. The tail wind case has a +5.8 m/s (+19 ft/sec) di^fference in 
airspeed which causes - 14 sec time and -1.8 kg (-4 lb) fuel burn differences. 
These differences could possibly be lessened with airspeed control added. 


Table 5. Comparison of Optimization and Trajectory Generation 
Program Results. Descending from 10668 m (35000 ft) to 1.8 m 


(6 ft). 


Type 

Fuel Bum - 

I - 

Tlme-to- 

-go - 8 

Altitude 

m 

■ (ft) 

Airspeed 

o/s 

"(ft/a) 

Range- 

to-go 

Wind 

OPTIM 

TRAGEW 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

Jo Wind 

237 

(523) 

237 

(523) 

0 

8 

1.8 

(6) 

-3.1 

(-10) 

124 

(407) 

126 

(413) 

0 

-49 

(-157) 

Head 

230 

(506) 

234 

(516) 

0 

7 

1.8 

(6) 

-5.2 

(-17) 

124 

(407) 

121 

(396) 

0 

-96 

(-314) 

Tall 

243 

(535) 

241 

(531) 

0 

-14 

1.8 

(6) 

-4.3 

(-14) 

124 

(407) 

130 

(426) 

0 

-85 

(-280) 


Despite the differences between OPTIM and TRAGEN results, it is believed 
that the match between them is very good. That is, the results provided by 
OPTIM (fuel burn, time expired, trajectory followed) can be concluded as 
being accurate, and sensitivity runs based on using both programs will also 
produce accurate conclusions. A better match between the two programs can 
be obtained by minor adjustments to both programs. 


Sensitivity Analysis 


In conducting the sensitivity analysis of the optimum profiles, there 
are two items that are of particular interest: 

1) . The establishment of how the characteristics of optimum profiles 
change with variations in initial weight, range-to-go, wind pro- 
files, and other variables affecting the aircraft performance. 
Knowing these variations affects how an on-board algorithm should 
be constructed to account for measured changes in the flight 
conditions . 
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2) . The analysis of the effect that errors In the assumed flight con- 
ditions (i.e., different initial weight, different wind profile) 
have on the cost of flying a particular profile. For example, 
if assuming the wrong initial weight of the aircraft has little 
effect on the overall cost, then the algorithm can be constructed 
so that initial weight is assumed to be a nominal value. 

The purpose of this section is to begin to address these items for wind, 
initial weight, and range variations. Both OPTIM and TRAGEN are required 
for this purpose. Other variables that may have an effect are temperature, 
lift coefficient, drag coefficient, and thrust variations as well as 
measurement errors of altitude, flight path angle, airspeed (or Mach 
number), and range-to-go. These variations should be addressed at a future 
time. 


It is assumed for climb that Mach number and flight path angle as 
functions of altitude specify the optimum profile to a given cruise altitude. 
Thus, by presenting plots of these variables as functions of altitude, the 
effect of the variations to the nominal conditions is directly seen. 

Figure 10 shows the variations in the climb profile as the total range 
is varied, in steps of 25 n.mi. from 75 n.mi. to 275 n.mi. As can be seen, 
there is no variation below 3048 m (10000 ft) . Above 3048 m (10000 ft) , the 
flight path angle is initially higher and the climb Mach number is initially 
lower as the total range is increased. But it is seen that these curves 
are essentially parallel, and they merge into a common curve for range ex- 
ceeding 175 n.mi. Thus, it would be relatively easy to model these curves 
as a polynominal function of altitude with parameter variation due to change 
in range. 

Figure 11 shows the variation in the profile as initial weight is 
varied. Again, there appears to be a parallel offset with a greater effect 
on flight path angle. These is no effect on Mach number below 3048- m 
(10000 ft).) Also, again it can be seen that this variation would be easy 
to account for in specifying an on-board reference trajectory. 
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Mach Number Flight Path Angle 




Mach Number Flight Path Angle 



The biggest uncertainty that would be experienced during the climb 
and descent would be the variations in the longitudinal wind. To study 
this effect, the two wind profiles of Fig. 12 were used as program inputs. 
The Denver wind is the same as is shown in Fig. 9, but Fig. 12 also has 
the wind headings indicated at discrete points. The "triangle" wind is 
a constant 270° heading wind varying linearly in magnitude with altitude 
at a rate of 1 knot per 305 m (1000 ft) . These winds profiles were assumed 
to act as both tail and head winds by setting the aircraft heading to be 
90° or 270°. 

Figure 13 shows the effect of each of these wind profiles on the 
optimum climb trajectory compared to the nominal profile having no wind. 

As can be seen, the Mach number variation with altitude is essentially 
a parallel offset. The effect on flight path angle is mainly seen going 
from zero to 3048 m (10000 ft) where a fan out of +0.3° is seen. Above 
3048 m (10000 ft), the variation appears more as a steady offset. Thus, 
again it appears that accounting for various wind profiles can be done 
in the on-board system in a relatively simple way. This conclusion is 
discussed more in the next chapter. 

Figure 14 shows the effect of each of the wind profiles on the optimum 
descent from 10058 m (3000 ft). There, the plots show altitude and Mach 
number as functions of range-to-go to the landing point. From the first 
plot, it is seen that wind causes a + 10 n.mi. range variation in where 
the optimum descent should begin. This is predictable based on knowing 
how far the air mass moves due to the wind during the descent period. 

Also, there is up to a + 0.05° variation in the optimum Mach number as 
a function of range-to-go. From the earlier results presented, it was 
shown that if the altitude-range profile was followed closely down to 
3048 m (10000 ft), the Mach number would also follow the nominal profile. 
Thus, from an implementation point-of-view, the Mach variations seen 
above 3048 m (10000 ft) in Fig. 14b are not of concern. Below 3048 m 
(10000 ft), only one profile needs to be stored. 


41 



Altitude 



Figure 12. Wind Profiles Used to Study Sensitivity Effects on Optimum 

Profiles. 
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The results seen in Figs, 10''14 are limited in scope, in that all possible 
variations are not addressed. For example, the effects of variations in 
cruise weight and cruise altitude for the descent trajectories should be 
explored. However, the one thing that can be concluded is that variations 
of total range, initial weight, and wind profiles have predictable effects 
that can be easily used to modify the characterization of the nominal 
optimum trajectory. Thus, the process of computing the optimum climb-des- 
cent profiles does not have to take place on-board the aircraft. Instead, 
a nominal profile plus modifications to account for off-nominal parameter 
changes can be pre-computed and stored in the aircraft flight management 
system. This provides a very simple on-board method of computing the 
optimum profile. 

The second aspect of the sensitivity analysis was to study the effect 
of flying a non-optimal profile. This effort involved using TRAGEN to 
simulate following a given incorrect optimum profile when a different 
profile should have been used. Two error conditions were investigated; 

a) Following an optimum profile specified for the incorrect 
initial weight. 

b) Following an optimum profile specified for the incorrect wind 
profile. 

Only the climb phase was investigated. 

For initial weight errors, four runs were made as shovm in Table 6. As 
can be expected, the increase in assumed initial weight causes a decrease 
in the optimum cruise altitude. Thus, if the aircraft is lighter than the 
assumed amount, it will climb faster to the lower altitude. This results 
in both time and fuel reduction from what was predicted. Cases la and Id 
in Table 6 were of this nature, and they show a 13-14% reduction in fuel 
and a 14-16% reduction in climb time for an initial error of -6804 kg 
(-15000 lb) in initial weight. 
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Table 6. Initial Weight Error Cases Considered - Climb Profile 


Run 

Initial Weight 
- kg (lb) 

Assumed Optimum Profile i 
Weight - kg (lb) 

Change in 
Climb Fuel 
kg (lb)-(%) 

Change in 
Climb Time 
Sec - (%) 

la 

61236 

68040 

-774 



(135000) 

(150000) 

(-604) (-131%) 

-159 (-14.4%) 

b 

61236 

54432 

+329 



(135000) 

(120000) 

(+726) (+18.7%) 

223 (+22.9%) 


68040 

61236 

371 



(150000) 

(135000) 

(+817) (+19.2%) 

247 (+23.3%) 

d 

54432 

61236 

269 



(120000) 

(135000) 

(-592) (-13.9%) 

-167 (-15.8%) 


On the other hand, if the initial weight is assumed too small, the 
heavier aircraft will attempt to climb to a higher than nominal altitude. 
Cases lb and Ic of Table 6 show this condition for an initial weight 
error of +6804 kg (+15000 lb). The result was a 19% increase in required 
climb fuel and a 23% increase in required climb time. It was also seen 
that the aircraft did not achieve the required cruise airspeeds when the 
incorrect optimum input altitude was reached. 

Based on the results of Table 6, it is seen that the initial weight is 
an Important parameter to be entered into the computations. It is better 
to assume too large an initial weight than vice versa. 

For the climb wind profile errors, four more cases were studied where 
optimum climb profiles based on the triangle wind profile of Fig. 12 were 
used. The results are shown in Table 7. As can be seen, if the tall wind 
is greater than assumed (Cases 2a and 2d), more time and possibly more fuel 
is required to achieve cruise conditions. Likewise, if the head wind is 
greater than assumed (Cases 2b and 2c), less time and possibly less fuel 
is required, although the results are mixed. It is seen that neglecting the 
wind profile can vary fuel cost up to + 1.5%. However, this effect is not 
nearly as significant as assuming the incorrect value of aircraft weight. 
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Table 7. Wind Profile Error Cases Considered - Climb Profile 


Run 

Actual Wind 
Profile 

Assumed Profile 
for Computing 

Change in 
Climb Fuel 

Change in 
Climb Time 
Sec - (%) 

Optimum 

kg (lb) - (%) 

2a 

No Wind 

Head 

-1 (-2) (-.05%) 

+7 (+0.7%) 

b 

No Wind 

Tail 

+11 (+24) (+0.6%) 

-3 (-0.3%) 

c 

Head 

No Wind 

-27 (-60) (-1.4%) 

-15 (-1.4%) 

d 

Tail 

No Wind 

+29 (+64) (+1.5%) 

+18 (+1.7%) 


These type of sensitivity studies are important because they produce 
information necessary for the implementation process. Computing the 
optimum profile with the wrong aircraft or environment models can cause 
a large percentage of the expected gain to be obtained from the flight 
management system to not be realized. 

Many more sensitivity cases than those described above need to be 
obtained for determining the sensor and measurement processing requirements 
associated with Implementation of an optimum vertical flight management 
system. However, with the availability of the OPTIM and TRAGEN programs, 
the user is in a position to obtain these results. 
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V 


MECHANIZATION OPTIONS 

This chapter briefly describes two methods by which the optimum 
vertical profiles can be generated and used for automatic longitudinal 
steering on-board the aircraft. These are: 

1) Computing the nominal optimum profile off-line and storing it 
on-board the aircraft. Empirical equations are used to adjust 
this nominal profile to account for perceived changes in initial 
weight, range, wind profile, etc. 

2) Continually recomputing the nominal optimum profile on-board 
the aircraft. Here, changes in wind profile and range are 
continually being estimated, and their effects are directly 
entered into the nominal trajectory computation. 

Each method is now discussed. 

Off-line Optimum Profile Generation 

This method would be based on sensitivity analysis results such as 
those presented in the previous chapter. A nominal optimum profile is 
computed off-line using the OPTIM program described earlier. Also, several 
other profiles are computed, where variations are made (one at a time) in 
the various key parameters that govern the shape of the profile. Parameters 
that would be varied include initial weight, wind profile, temperature 
profile, and range- to-destination. 

A means would be derived for describing the nominal profile as a 
function of a suitable independent variable. For example, it is convenient 
to specify airspeed (or Mach number) and flight path angle as functions of 
altitude for climb and altitude and airspeed as functions of range-to-go for 
descent. Also, empirical functions would be derived for computing perturba- 
tions to the nominal profile for measured changes to nominal values of the 
key parameters. 
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The function describing the nominal trajectory and its perturbations 
are stored in the on-board flight management system's computer. Then, as 
the chosen independent variables change with the position along the flight 
path, the sum of the nominal and perturbation functions are used to compute 
and command the desired flight profile. 

An example of this mechanization for the climb profile is depicted in 

the block diagram of Fig. 15. Here, the independent variable is altitude 

h. Its measurement is used to compute the nominal values of airspeed 

V^^(h) and flight path angle The Nominal Trajectory computations 

are also used to compute the nominal wind profile V^(h) and the assumed 

aircraft weight W (h) . Other measurements are used to estimate or measure 
^ ^ 

the prevailing values of wind and weight W. If range-to-destination 
r^ is changed, this is also entered. These variables are compared to the 
nominal values in the Empirical Perturbation Functions computations. 

These computations produce perturbations 6V^(h) and 6y^(h) to the nominal 
commands. The perturbations and nominal values are added to obtain the 
total commands V^(h) and Yj,(h). 



Figure 15. Block Diagram Depicting On-Board Climb Steering Based on 
Computing the Nominal Optimum Profile Off-line. 
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The total commands are compared with measured or estimated values of 
airspeed and flight path angle (V, Y). The differences ((SV,6 y) used 
to implement a linear control law based on perturbation analysis. This 
control law produces either a commanded angle-of-attack, pitch angle, or 
elevator deflection to produce the desired longitudinal response character- 
istics. The gains in the control law could be functions of altitude if 
desirable. Maximum thrust would be set during climb. 

A similar mechanization would be used for the descent portion of 
flight but with the independent variable changed to range-to-go. For 
descent, the thrust would be set at idle value. 

For cruise, the thrust is set to balance the drag at the best cruise 
airspeed. This airspeed is chosen to minimize the cost per unit of dis- 
tance traveled. If a constant cruise altitude is required because of ATC 
constraints, then the elevator is trimmed so that lift balances the air- 
craft weight. 


On-Board Optimum Profile Generation 

This method would be based on the algorithms used to generate the 
optimum profile, as in the OPTIM program. The algorithms would be exercised 
every few minutes as background equations in the on-board flight management 
computer. The sequence of background computations would be similar to the 
seven steps described in Chapter III. As new estimates of the wind profile, 
aircraft weight, or other key variable were obtained, a new optimum pro- 
file would be generated based on the latest information. There would be 
no need for perturbation functions . This method would require a much larger 
computational capability then the first method described above. Information 
required by the system would include : 

1) minimization algorithm; 

2) functions or tables describing: 

a) standard atmosphere. 
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b) aircraft lift characteristics, 

c) aircraft drag characteristics, 

d) engine thrust characteristics (maximum and idle) , 

e) fuel burn characteristics; 

3) function interpolation schemes; and 

4) empirical equations for estimating cruise weight and 
landing weight. 

An example of this mechanization for the climb portion of flight is 
shown in Fig. 16. Note that a new nominal trajectory is computed at 
reqular intervals. After the commanded values of airspeed and flight path 

angle (V^(h), are computed, the mechanization is identical to that 

shown in Fig. 15. 

This type of mechanization has been developed for a DC-10 and exercised 
in a cockpit simulator, as reported in Ref. 25. For short range flights 

(200 n.mi.), this mechanization was reported to save over 13% of the fuel 
requirements . 


A 



Figure 16. Block Diagram Depicting On-Board Climb Steering Based on 
Computing the Nominal Optimum Profile in the Background. 
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Instrumentation required to support these two flight management systems 
includes: 


a) 

Vertical accelerometer ^ 


b) 

Airspeed measurement system (Mach and/or true 

airspeed) ^ 

c) 

Baro-altimeter , 


d) 

DME (or INS) so distance-to-go and wind speed 

can be obtained. 

e) 

Free stream temperature measurement, and 


f) 

Vertical gyro to measure pitch angle. 



A complementary filter can be used to blend vertical accelerometer and baro- 
altitude measurements to obtain altitude rate. INS and airspeed measurements 
can be combined to obtain wind estimates. Altitude rate and airspeed can 
be combined to compute flight path angle. Pitch angle and flight path angle 
can be differenced to obtain angle-of-attack. Angle-of-attack can be used 
with temperature, airspeed, and altitude to compute lift and aircraft 
weight. 

Also, this system must be coupled to measurement of the throttle and 
elevator deflections. 


Method Comparison 


The two methods described above cannot be directly compared without 
actual mechanization on the same type of computer and then tested under 
realistic operating conditions. However, some qualitative comparisons 
can be made before that point. Advantages to each method are: 


Off-line Method 


1. Simpler computations on-board. Does not require table interpola- 
tions, iterative optimization procedures, or foreground/background 
time sharing. 
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2 . 


Smaller memory requirements. Few subroutines are necessary. 
Tabular data minimized. 


3. Easy to incorporate simplifications • 

4. Less expensive. 

5. No concern over converging to correct answer. 

On-board Method 

1. More accurate optimum profile generated. 

2. Need for flight engineer to enter information before flight 
is minimized. 

3. More adaptable to changing flight conditions (i.e. re-routing 
due to weather constraints) . 

4. Not restricted to variations in a limited number of key 
parameters . 

Thus, it is seen that both methods have advantages. It is recommended that 
both of these methods be investigated further by mechanizing them on an 
airborne system and then testing them in flight (or in a cockpit simulator 
using a realistic operating environment) so that further evaluation may 
be made. 
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VI 


SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS 


Summary and Conclusions 


All evidence indicates that we are rapidly consuming our supply of 
hydrocarbon fuels. This makes it mandatory that systems be developed that 
allow conservation of the remaining supplies until alternate fuels can be 
developed. This is especially true of aircraft flight which currently 
has no alternatives. This report addressed the development of on-board 
algorithms for vertical steering of the aircraft to minimize fuel con- 
sumption and cost. 

Chapter II presents a summary of the research work conducted in the 
past twenty years to develop algorithms for fuel minimization and related 
topics. The trend has been to use the energy state approximation for 
obtaining a simple mathematical model of longitudinal aircraft dynamics. 

This simplification has facilitated development of algorithms that are 
relatively easy to solve so that they can be mechanized on-board the air- 
craft for in-flight solution. Furthermore, based on several researchers' 
results, fuel savings of from 5% to over 25% can be obtained by using these 
algorithms. (See Refs. 3, 4, 25, and 26). 

Chapter III derives an algorithm for computing the optimum vertical 
profile using range as the Independent variable. Both fuel and time are 
penalized, and the longitudinal wind effects are taken into account. The 
Hamiltonian is constant for this mechanization, and is equal to the mini- 
mum cruise cost per unit distance traveled. To obtain optimum climb and 
descent profiles Involves minimizing a single function at discrete points 
along the trajectory by proper choice of thrust and airspeed. This algorithm 
proved to be a dual to the one derived by Erzberger where energy was used 
as the independent variable. A computer program (called OPTIM, described in 
Appendix A) was used to obtain optimum vertical profiles' for typical 727 
aircraft flights based on these algorithms. 
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In Oiapter IV, the accuracy of the vertical profiles obtained from 
OPTIM were examined by using a more complete longitudinal model of the 
aircraft. This model was incorporated into a computer program called 
TRA.GEN (described in Appendix B) which steers the aircraft to follow the 
input reference trajectory. Results showed that the OPTIM reference tra- 
jectories were both flyable and accurate in terms of the fuel burned and 
time expended. 

Chapter IV presented the changes in characteristics of optimum 727 
reference trajectories due to changes in range, initial weight, and wind 
profile. It was concluded that these changes were simple modifications 
to the nominal reference profiles; they could be used to compute perturba- 
tions to a nominal profile on-board without recomputing the entire reference 
profile. 

Chapter IV also examined the effects of errors in the estimated values 
of initial weight and the wind profile on the performance obtained during 
climb. A 9% increase in initial weight (6804 kg (15000 lb)) can cause a 
23% Increase in time and a 19% Increase in fuel required to achieve ‘the 
desired cruise conditions. Wind errors have a smaller effect. These sen- 
sitivity studies are useful for specifying how accurately various para- 
meters which affect the flight performance need to be measured. 

Chapter V presented two alternate ways that optimum vertical profiles 
could be computed and mechanized in flight . One method is based on com- 
puting the reference profile and perturbation adjustments off-line, so that 
minimum computations are required on-board. The other method is based on 
continually computing the optimum reference profile on-board the aircraft 
so that measured changes to the nominal conditions can be directly taken 
into account. Both methods have qualitative advantages, and a choice 
between the two would have to be based on further mechanization studies. 
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Recommendations 


Based on the experience obtained by review of previous work, the use 
of the OPTIM program, the development and use of the TRAGEN program, and 
conversations held with many government and industry personnel associated 
with research and design in the air transportation field, the following 
recommendations are made concerning future work. 

OPTIM and TRAGEN extensions . The following extensions to these pro- 
grams would provide the capability to obtain more realistic reference 
trajectories, correct known minor errors, and add capabilities which would 
be useful for obtaining improved solutions for other phases of flight: 

a) Optimum climb to a fixed altitude . Currently, the aircraft is 
simulated to climb to the optimum altitude without considering 
ATC constraints. Also, OPTIM computes cruise costs based on the 
assumption that the aircraft will climb during cruise as fuel 

is burned. This added option would constrain the profile to 
climb to a fixed altitude and then hold it during cruise. 

b) Flight path angle smoothing . Currently, because OPTIM computes 
climbs and descents at discrete changes in specific energy, 
jumps occur in the flight path angle producing a rough reference 
trajectory. Airspeed change constraints and post-run smoothing 
should be tested to remove or minimize the effect of this problem. 

c) Optimum profile with multiple cruise minimums . If there is a 
high head wind shear condition, there will exist more than one 
best cruise altitude depending on range. Currently, OPTIM will 
only consider the lower altitude. Thus, for longer ranges, with 
certain wind profiles, better flight conditions than what OPTIM 
will produce will exist. The change would add logic to consider 
all altitude minimums in deciding the cruise condition. 
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d) Minimum fuel with fixed time of flight (4D requirement) . OPTIM 
produces a profile which minimizes the cost function 

J = 

o 

where and are input costs of fuel ($/lb) and time ($/hr) . 
OPTIM can be set up to Iterate on (with held fixed) so 
that the desired flight time would be produced for 4D guidance. 

e) Optimum cruise - step change in altitude . Currently, OPTIM is 
based on optimum cruise climbs. This added option would be an 
extension to (a) above where cruise altitudes are fixed by ATC 
constraints. This option would determine where to change altitudes 
in a step fashion as fuel is burned off during cruise. 

f ) Optimum cruise change in altitude due to range dependent wind 
profile change. For certain aircraft on certain routes in 

the future, the ATC requirement for fixed altitude during cruise 
will be removed. However, for long flights, the wind profile 
encountered during cruise is continually changing, and it will 
be desirable to adjust cruise conditions continually to minimize 
costs. This problem may require transient engine thrust and 
fuel flow models as a function of variable throttle. Also re- 
quired are typical models of the changing enroute wind profile. 
OPTIM could be used to generate optimum cruise tables for vari- 
able wind and cruise weight. How to climb and descend as a 
function of range would be determined. 

g) Follow arbitrary profiles. Currently, TRAGEN is designed so that 
the aircraft tracks a sequence of input points which come from 
OPTIM. Equations could be added so the aircraft could also track 
constant indicated airspeed, constant Mach number, profile des- 
cents, and other possibilities of interest. The result would be 
the ability to measure the cost (in time and fuel) to fly to 
follow these alternate reference profiles. 


I 


(CfW 


+ C 


t)dt 
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h) Maintain airspeed during descent. Currently, TRAGEN is programmed 
so that angle-of-attack is changed to track inertial flight path 
angle (i.e., altitude as a function of range) during descent. 

Thrust is set to idle value. Below 3048 m (10000 ft), the airspeed 
drops below 250 kts causing excess time and fuel burn. An option 
is to add throttle control to maintain appropriate airspeed. 

Another option is to add airspeed error to the control law as is 
done for climb control. 

i) Conversion of OPTIM and TRAGEN to model the Boeing 737. Currently 
OPTIM and TRAGEN are based on mathematical models of the 727 
aircraft. Because future research will be focused on the NASA 
Langley 737 TCV aircraft, it is appropriate to develop an alter- 
nate version of these programs based on 737 models. The sub- 
routines and data tables associated with aerodynamic and propul- 
sion characteristics would have to be changed. Also, empirical 
equations which exist within OPTIM associated with estimating 
initial cruise weight and final landing weight would have to be 
changed. This task would also include a written explanation of the 
steps involved in making the conversion. 


Develop and test optimum vertical guidance methods. The purpose of 
this task would be to derive the specific requirements to design and test 
an airborne system for implementing optimum vertical guidance techniques 
on either a 737 cockpit simulator or the TCV aircraft. Both methods des- 
cribed in Chapter V should be considered. Subtasks include: 

a) Generate perturbation guidance functions (Method I*) This task 
would parameterize the optimum profiles for the 737 (offline, as 
determined by converted OPTIM) to account for measured or input 
changes in range-to-destination, takeoff weight, wind profile, or 
temperature profile. This would involve describing the climb 
and descent profiles as functions or polynomials of altitude, 
range, etc. These parameterized profiles would be used to 
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generate the total (nominal plus perturbation) reference trajec- 
tories for climb and descent steering. Control gains would be 
scheduled as necessary to give appropriate steering response. 

b) Test perturbation guidance on TRAGEN. This would involve testing 
the perturbation steering commands in the presence of wind, weight, 
temperature, range, and other variations on TRAGEN. The results, 
in terms of time and fuel expenditures, should match the results 
predicted by OPTIM. If not, the parameters in the profile equa- 
tions would be adjusted. The result should be a set of near- 
optimum steering equations for the 737 that are simple to imple- 
ment yet complex enough to account for most perceived perturbations. 

c) Generate on-board optimum profile guidance computations (Method 2) . 
This task would specify and code the tables, logic, and equations 
required for direct on-board computation of the optimum vertical 
profile. Essentially, this would duplicate the OPTIM calculations 
on-board the 737. This code would be validated by using the 
TRAGEN program. 

d) Determine airborne computer mechanization requirements . The com- 
puter where the optimum guidance algorithms would reside would be 
chosen. The interface between the computer and the sensors and 
autopilot (or flight director) would be defined. Also, the inter- 
face between the vertical profile steering algorithms and other 
avionics functions would be defined. Block diagrams and flow 
charts would be developed of the required equations and logic 

in sufficient detail so that coding could take place. 

e) Compare performance and mechanization requirements. Based on the 
results of subtasks (a) - (d) above, the performance and require- 
ments for mechanizing the perturbation guidance algorithm (Method 1) 
would be compared with the performance and mechanization require- 
ments of the optimum profile guidance algorithm (Method 2). The 
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comparison would include computation speed and memory require- 
ments. Based on these results, a decision could be made con- 
cerning further mechanization on the TCV aircraft or 737 cock- 
pit simulator. 

f) Mechanize and test the preferred algorithm. Based on subtask 

(e), the chosen algorithm should be mechanized. A test plan should 
be developed, an extensive test should be made, and the results 
should be publicized. In this way, all avionics firms, airlines, 
and personnel associated with air transportation would have 
access to a near-optimum method for saving fuel and reducing 
flight costs. 
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APPENDIX A 


TRAJECTORY OPTIMIZATION USING THE ENERGY STATE METHOD 

The purpose of this Appendix is to summarize briefly the theore- 
tical background used for optimization of trajectories. These principles 
form the basis of the numerical process used for computing the optimum 
vertical profile of a jet aircraft. More details are given in Refs. 5, 

24, and 30. Reference 5 derives the principles upon which trajectory 
optimization is based. In Refs. 24 and 30, Erzberger and Lee apply these 
principles using the energy state approximation to obtain a practical, 
efficient means of generating the optimum vertical profile. 

OPTIM is an extension of the original computer code developed by 
Erzberger and Lee that encompasses these principles and is based on their 
methods. The latter part of this appendix describes how OPTIM is organized. 

In the following sections, the theory of trajectory optimization is 
first presented. Then, the application of this theory to minimizing the 
direct operating cost (DOC) of an aircraft traveling over a fixed range 
is outlined. This is followed by a discussion of the details of going 
from theoretical expressions to a practical computer code. The theore- 
tical points are presented without proof, for conciseness. The reader 
wanting more detail should review the references. 

Theoretical Principles 

In Ref. 5., a description is given of the requirements for solving 
an optimization problem involving a continuous dynamic system with no 
terminal constraints but with fixed terminal time. This description is 
repeated here because it presents the basic principles which extend to 
the aircraft profile optimization problem- 
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A system (the aircraft) is governed by the nonlinear differential 
equations 

X = f(x,u,t) ; given; (A.l) 

^o 1 ^ 1 ^f ’ 

where x is the n-dimensional state vector and u is the m-dimensional con- 
trol vector. The cost function which is to be minimized is of the form 


<{>(x(t^) ,t^) + J' ^L(x,u,t) dt. 


(A.2) 


Here, ((> is the terminal cost function, and L is the cost per unit time along 
the trajectory. The problem is to find the sequence of controls u(t) that 
minimize J. 


First, the system equations are adjoined to J with the multiplier 
vector X(t) : 

/ ^f T • 

{L(x,u,t)+X (t) {f (x,u,t)-x}} dt. (A. 3) 

*"o 

Then the Hamiltonian function is defined as 

H(x,u,t) = L(x,u,t) + x'^Ct)! (x,u,t) . (A. A) 

Equation (A. 3) is integrated by parts to yield 

J = 4>(x(t^),t^) - X^(t^) x(t^) + X^(t^) x(t^) (A. 5) 

+ ^ ^ {H(x,u,t) + x'^(t) x(t)} dt. 
t 

o 

Next, the change in J due to variations in u(t) and x(t) is con- 
sidered for fixed t and t,: 

o f 
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6J. 


J 5x I + (X^6x) 

/ ) t=t . ■ o 


+ |f 6u dt. 


(A.6) 


The elements of A(t) are chosen to cause the coefficients of 6x in Eq. 
(A.6) to vanish under the integral and at t^: 

?T 9H 9L ,T 9f ,T. , 9d) 


(A. 7) 


Equations (A. 7) are called the co-state equations. Then, Eq . (A.6) becomes 


6J = X 6x(t^) + 


6u dt. 


(A.8) 


For J to be minimum, 6J must be zero for arbitrary u(t); this Implies that 
for no bounds on u, 


■^=0 , t<t<t. 

3u ’ o — — f 

on the optimum path. If the control variables are constrained as 


(A.9) 


C(u,t) £ 0, (A. 10 : 

then for u(t) to be minimizing, we must have 6J £ 0 for all admissible u(t) 
This implies, from Eq. (A.8) that 

1^ 6u ^ 6 h > 0, (A. 11' 

du = — 


for all t and all admissible 6u(t). In other words, H must be minimized 
over the set of all possible u; this is known as the minimum principle [6]. 


In summary, to solve for u(t) that minimizes J, the differential 
equations (A.l) and (A. 7) must be solved simultaneously, where u(t) is 
determined from Eqs. (A.9) or (A. 11). The boundary conditions on the state 
X at t^ and X at t^ are specified, resulting in a two-point boundary-value 
problem. 
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If L and f are not explicit functions of t, then 


H 


M 4. M * 4. M 

3t 8x ^ 9u 


u + 



8H 

3t 


4. * 4. 




(A.12) 


Each element of Eq. (A.12) is zero on the optimum trajectory, from which 
we can conclude that H is constant on the optimum trajectory. This latter 
point is used in the analysis presented in Chapter III. 


Application to Aircraft Profile Optimization Using the Energy 

State Approximation 


Here we are concerned with applying the above theory to the pro- 
blem of choosing the thrust and airspeed values to control the aircraft 
vertical profile in going from one point to another. The cost function 
J is the direct operating cost (DOC) which is the sum of fuel and time 
costs. This is, in integral form, 

J = w + Cj.)dt f (A. 13) 

where C^ is the cost of fuel ($/lb), w is fuel flow rate (Ib/hr), C^ is the 
cost of time ($/hr), C^ is the direct operating cost, and t^ is the 
time to fly the specified distance traveled d^. It is also assumed that the 
typical vertical profile is as shown in Fig. A.l - that is, it contains 
climb, cruise, and descent portions which have the constraint that 


d 

up 



(A. 14) 


where 


d = x(t .) = the distance traveled from the start point to 
up ci 

where the cruise segment begins (at time t = t ,). 


d, = d^ - x(t ,) = the distance traveled from the end of 
dn f ' cf 

cruise (at time t = t^^) to where the descent segment ends. 
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Thus, the cost function (Eq. (A. 13)) can be rewritten as 

j =y'""c^dt+(d^-d^p-d^^),j; + c^dt, (A.15) 

0 

= c^dt + (df- x(t^.)- [d^ - + y^f . c^dt. 

*0 1 4: 

cf 

where is the cost per unit distance while in cruise. 

Simplified point-mass equations of longitudinal motion of the air- 
craft are 


V = (T-D)/m - g sin Y , (A. 16 ) 

h = V sin Y , 

X = V cos Y + V , 
w 

where the flight path angle (Y) dynamics and weight loss due to fuel burn 
are neglected. Here, 

V = airspeed, 
h = altitude, 

= longitudinal component of wind speed, 
m = aircraft mass, 

T = thrust, 

D = drag. 
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Also, it is assumed that lift L = mg cos Y* The effect of weight loss is 
accounted for by continuously updating weight without adding another state 
variable. 

Specific energy E is defined as 

E = h + V^/2g, (A.17) 

which is the sum of potential and kinetic energy per unit mass. Its 
time derivative is found to be 

E = V(T-D)/mg. (A. 18) 

The energy state approximation is based on the assumption that potential 
and kinetic energy can be Interchanged instantaneously. In this approxima- 
tion, the energy state variable replaces altitude and airspeed state 

• • 

variables [12]. Thus, Eq. (A.17) can be used in place of V and h in Eq. 

(A. 16). 

It is assumed that the aircraft specific energy Increases monotonically 
during climb and decreases monotonically during descent. This assumption 
is used in the development to change the independent variable in Eq. (A. 15) 
from time to energy. This uses the transformation 

dt = . (A. 19) 

E 

It is mathematically convenient to evaluate the last integral in Eq. 

(A. 15) backwards in time so that the energy state is monotonically increas- 
ing during its evaluation. This means that the running distance (range) 
variable during the descent can be measured backwards from the end point. 
Thus we can think of range measured in two ways as shown in Fig. A. 2. 
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CRUISE 


CLIMB, 


DESCENT 




X , ( t - ) - 

dn cf 


Figure A. 2. Measurement of Range from the Origin or to the Destination. 


In this sketch. 


X (t) = 

up 

*up<'cl)= 


range measured on the way up In forward time t, 

value of X when initial cruise is reached, 
up 

range measured on the way up in backward time T, 

value of Xj when final cruise is reached (t ^“It ,-t-|). 

dn cf ' cf f ' 


Also, we define the variable x to be range traveled during climb and 
descent. The distance traveled during cruise is then contrained to be 
(d^ - x) . We can then see that an incremental change dx in the range 
variable x is equivalent to Incremental changes in both x^^ and 
That is 

dx = d(x + X , ) . ( 

up dn“^ 

From this discussion, the second of Eqs. (A. 15) can be written as 

J -/‘clVt 4. (d, - 4 fcl IcJdT . 


(A. 20) 


(A. 21) 


We use Eq. (A. 19) and the transformation 


(A. 22) 
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to rewrite Eq. (A. 15) as 





+ 



(A. 23) 


Here, E^, are the values of specific energy evaluated 

at time t equal to 0 and t^^ and time t evaluated at t^^ and respec- 
tively. 


Note from Eq. (A. 23) that the range variable x only appreas as the 
sum of cli m b and descent distances + ^dn^ * state equation 

for this system of equations can be written as 


dx 

dE 


d(x + Xj ) 
up dn 

dE 


/(V + V ) \ / 

\ E / E>0 ' 


(V , + V , 
dn wdn 




E<0 
(A. 24) 


Here, and are the longitudinal components of the wind speed 

for climb and descent. Then, analogous to Eq. (A. 4), the Hamiltonian is 


H = 


( t ). ^ iU ^(™). ' 

E>0 ' ' E<0 E>0 ' E<0 


(A. 25) 


This can be divided as 


H = 


• « 



C , + X (V +v ) 
d up wup 

+ 

E>0 


• 

E 

|£| 


(A. 26) 


J E<0 

Now, analogous to Eq. (A. 7), the costate equation for X can be written as 


70 



^ ^ _ IM = _ iH 

9[E 8x d(x +x, ) 

up dn 


(A. 27) 


and from Eqs. (A. 7) and (A. 23), this costate has the final value 


X(E ,) = X(E .) = Tr— Z T 

ci cf 3(x +x, ) 

up dn 


9([d.-x -X, ]ij^) 

f up °n ^ , 

3 (x +x ) ^ 

''^up ^dn-* 


(A. 28) 


where ip is the cruise cost per unit distance. 


Note, this problem could be placed in a slightly more conventional form 
by dividing it into two problems - one for climb and one-half of the 
cruise distance and the other for decent and the other half of the cruise 
distance. Then Eqs. (A. 27) and (A. 28) would be replaced by 


3e 3x 


A(E .) 

Cl 


3([d./2 - x ] ip(E .)) = - ip (E .), 

f up Cl Cl 

3x 

up 


for climb. For descend. 


3X ^ 3H 
3E 3x 

dn 


X(E^j) - 


3([dj/2 - 


3x 


dn 


= - ^ (E,f) 


(A.29) 


(A. 30) 


This allows splitting the Hamiltonian defined in Eq. (A. 26) and allows for 
4 j(Eci) / 4'(E^^). In fact, in the actual implementation E^^ E^^ because 

optimum cruise energy changes as fuel is burned off. The principal results 
are unchanged, however. 

Thus, from Eq. (A. 11), (A.29) and (A. 30) the trajectory optimization 
problem becomes 
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H 

up 


min 

V 

up 

TT 


up 



- ^(E .) 

^ cx' 



min 

^dn 


dn 



<l^(Ecf) 


/V + V 
I up wup 

' E 


) 


E>0 


» 



(A. 31) 


Thus, the optimization problem reduces to solving pointwise minimum 
values of the algebraic functions defined by Eq. (A. 31) during the climb 
and descent portions of the trajectory. 


Equations (A. 29) and (A. 30) are the transversality condition for the 

free final state problem (d + d, < d„) with terminal cost. Thus, the 

up dn f 

constant value of X for climb and descent is found to be the negative of 
the cost per unit distance for cruise. 

The cruise cost ip (= - X) is found by assuming that the aircraft is 
in static equilibrium during cruise (T = D) , and that 




min 

V 

c 


(V„ + 


V ) 
w 


(A. 32) 


In other words, for any cruise altitude, there is an optimum thrust and 

airspeed such that the cost per unit distance >I^(E^) is minimized. The 

optimum cruise cost as a function of cruise energy is typically of the 

shape shown in Fig. A. 3. Thus, there is also an optimum cruise energy 

E ^ where cruise cost ii(E ^) is minimized. If the range is long 
copt ^ copt 

enough so that there is sufficient range to reach optimum cruise energy 

E it should be done, and the cruise conditions should be set so that 

copt 
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•H 

p 



E 

copt 


Cruise Energy, E^ - m (ft) 

Figure A. 3. Optimum Cruise Cost as a Function of Cruise Energy 


For the case where there is no cruise segment (d^ = + d^^) , 

the cost function contains only integral terms. Then, the transversality 
condition yields X = - ijj(t^). That is, X would be the negative of 
'l^(t^), where is the optimum cost for cruising at the highest point 

reached on the climb trajectory. 


The optimum cruise energy is only specifically reached when 

there is range enough to climb to and descend from the optimum altitude/ 


airspeed values, where ^ minimum. For ranges less than this 

value, the maximum value of E that is reached is a free variable less than 

c 

the optimum value. Its choice is made to optimize the cost function of 


Eq. (A. 23). 


From Eqs. (A. 23) and (A. 25), one can write 



at E = E . This is 
c 


(A. 33) 
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(A. 34) 


H 

c 



= 0 




where d is the cruise distance, and H is the total value of H (H + H, ) 
c ’ c up dn"^ 

at the cruise point. Thus, Eq. (A. 34) can be used, along with other 
characteristics of \p and H, to determine the relationship between ip, E^, and 
d^. The Hamiltonian evaluated at E = E^ is the cost penalty to achieve a 
unit increase in cruise energy. For > 0, Eq. (A. 34) can be written as 

dc = -H^/0<l^/3E)g ^ g (A. 35) 

c 

Figure A. 4 shows the family of trajectories which have this characteristic. 

These occur at values of E below E ^ where 9ti^/9E < 0 (see Fig. A. 3). 

c copt 

That is, non-zero cruise segments occur at short ranges with cruise energies 
less than the optimum energy value for long range. 

For the case where H = 0, d is zero for d^/dE < 0. The distance 

c c 

d^ can be non-zero only at optimum cruise energy where 9i|j/9.E = 0. This 
family of trajectories is shown in Fig. A. 5. 


Thus, we have a situation where positive values of dictate one 

type of trajectory and zero values dictate another. In Ref. 24, it is 

shown that if the aircraft engine specific fuel consumption S__ is in- 

r L 

dependent of the thrust T (so that w = ST) , then the structure of the 

r L 

trajectories will be like Fig. A. 5 with no cruise segment occuring except 

at E (This implies that the Hamiltonian H is zero at the maximum energy 

copt ' r 

point). For this case, the optimum thrust setting for climb is and 

the optimum setting for descent is 

If the engine specific fuel consumption is dependent on thrust, and 

the thrust values are not constrained during climb or descent, it is 

shown in Ref. 2, that the Hamiltonian H is again zero at the cruise 

c 

energy, and again the trajectory structure is like those of Fig. A. 5. 
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If the is dependent on thrust, and constrained to the maximum 

r L 

value for climb and to the minimum idle value for descent, then the 
Hamiltonian is positive at cruise. This causes positive cruise segments 
according to Eq. (A. 35) at cruise energies below the optimum. For this 
case, the optimum trajectories will have shapes similar to Fig. A. 4. 

These trajectories are slightly less efficient than those of Fig. A. 5. 
because one less control is available for optimization. 

Some Mechanization Details of the Computer Program 

The remaining sections of this Appendix describe how the previous 
theoretical material has been utilized to construct an offline computer 
program for generating optimum vertical profiles for a model of the Boeing 
727-100 aircraft provided by the NASA Ames Research Center. This material 
is presented in an alternate way in Ref. 30, and the program is referred to 
here as OPTIM. 

By examining the specific fuel consumption data of the JT8D engine, 
it is determined that S_„ is dependent on thrust. Thus, for the 727, 
two types of short range profiles must be considered - those represented 
by Fig. A. 4 (Type 1 profile) when thrust is constrained and airspeed is 
the single control - and those represented by Fig. A. 5 (Type 2 profile) 

when both thrust and airspeed are used as controls. 

The solution to optimum climb and descent profiles is found by 
minimizing the Hamiltonian expressed in Eqs. (A. 31). The independent 
variable (energy) is stepped along in fixed increments (e.g., 150 m 
(500 ft)), and the Hamiltonian is minimized at each energy setting. 
Minimization occurs by finding the best values of air speed (V^p> 
and possibly thrust (F^p> so that the climb function and the des- 

cent function are individually minimized. 

To solve Eqs. (A. 31) requires knowing two more quantities: 

iKE^) - the cruise cost per unit distance. This comes from evaluating 
Eq. (A. 32) at the desired cruise altitude. • 
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E - the cruise energy. This is a function of the cruise altitude 
and the associated cruise airspeed obtained in Eq. (A. 32). 

Note that for the Type 2 profile at short ranges, there is no cruise 
segment. In this case, the maximum energy achieved at maximum altitude 
is referred to as the cruise energy E^. At that altitude, there still 
is defined a minimum cruise cost according to Eq. (A. 32). 


For the Type 1 trajectory of short range, there exists a non-zero 
cruise segment which is determined by use of Eq. (A. 35). To solve Eq. 

(A. 35) requires that the Hamiltonian defined by Eqs. (A. 31) be solved 
at the point of transition from climb- to-cruise. It also requires 
knowing the slope 3i^//3E of the cruise cost for a change in cruise energy at 
that point. 

Cruise Optimization 


The first step that must be taken to compute optimum trajectories 
is to derive the optimum cruise cost ip and its derivative 3if)/9E. This 
is done by computing what is referred to as the "cruise table". The 
parameters that affect this table are the assumed cruise weight, the 
wind profile, and the lift L, drag D, thrust T, and fuel flow rate w 
characteristics of the aircraft. The optimization process searches over 
the acceptable ranges of altitude and airspeed for a given weight. The 
results are collected in tabular form for a series of different assumed 
cruise weights. 


Again, the minimum cost of flight during cruise per unit distance 
for a fixed cruise weight is found by 




min 

V 


w + 
" V + V 



w 


(A.36) 
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This assumes that the aircraft is in static equilibrium during cruise, i.e., 


T cos a = D, 

L + T sin a = W, 


(A. 37) 


where the angle-of-attack a is found by solving these equations simul- 
taneously. The altitude is stepped in 305m (1000 ft) increments from sea 
level to ceiling altitude (where maximum thrust just balances drag). At 
altitudes below ceiling altitude, the airspeed - dependent drag curve 
crosses the maximum thrust curve at two points (Vj^ and V 2 ) as illustrated 
in Fig. A. 6. Thus, for each altitude level, the values of and V 2 
are determined, and then i(;(W^,E^) is minimized with respect to airspeed 
V between these two limits . Restrictions are that be greater than 
0.1 Mach and that V 2 be less than 0.89 Mach for structural reasons. 


After the cruise cost is minimized at each discrete altitude level, 

these numbers are stored in a table with altitude as the independent 

variable. Typical results are plotted in Fig, A. 7. Presented here are 

also the optimum cruise Mach number and the optimum thrust setting 

EPR After results are obtained in steps of 305 ni (1000 ft), the 

opt 

minimum cost point is found as a function of altitude. In the OPTlM 
program, the cruise table optimization results are obtained by using a 
Fibonacci search with 8 Fibonacci numbers. 


The cruise table results are obtained for cruise weights varying as 
dicated by the program input. Usually, the cruise weight is incremented 
in steps of 2268 kg (5000 lb) . Up to 10 values of cruise weight can be 
used. For each cruise weight, the optimal cruise altitude, cost, speed, 
power setting, fuel flow rate, and specific energy are computed. An ex- 
ample of optimum cruise cost as a function of cruise weight is shown in 
Fig. A. 8, 

Climb Optimization 

After the cruise tables are generated, the program procedes with 
obtaining the optimum climb trajectory. This requires guessing what the 
cruise weight will be, based on the takeoff weight. The guess is used 
to obtain a trial value for (or X) in the Hamiltonian from the cruise 
tables. The procedure to obtain this guess is based on an empirical for- 
mula which iterates until convergence is made. 
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Figure A. 8. Optimum Cruise Cost and Cruise Altitude as Functions of 
Cruise Weight for the 727 Flying into a Particular Head 
Wind. 


The climb optimization process starts by assuming 4' (E^) = (^copt^ ’ 

where '['(E is first obtained by setting the initial cruise weight W ^ 

copt ci 

equal to the takeoff weight (an input). The appropriate cruise tables 
are used to interpolate to find the corresponding value of associated 
with 1.5 4^. Then, an empirical formula of the form 


F = 0.11 (E , - E.) (1 + 4.7 C,/C.) W ./W . , 

up ci i t f Cl ref » 


(A. 38) 


is used to obtain an approximation to the fuel burned to reach E^. Here, 
E^ is the takeoff aircraft energy, is a reference weight (61690 kg 

(136000 lb) for the 727), and W^^ is the previous value of cruise weight. 
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Then, the cruise weight is updated at W . = W . - F . This process is 

cl Cl up 

repeated until the difference in consecutive estimates of F falls below 

up 

45 kg (100 lb) . 

When the cruise weight estimate is obtained, the feorresponding values 
of and obtained from the cruise tables. Then, the program 

is ready to generate points on the optimum climb trajectory. This is done 
by stepping along at discrete increments of specific energy and minimizing 
the Hamiltonian function 

C w + C - ) (V + V ) 

H (E) = — — : = — (A. 39) 

up ' E 

at each point. (This is the first of Eqs. (A. 31)). That is, the program 
starts with initial energy E^ = h^ + /2g. It steps the energy a fixed 

amount AE (say 150 m (500 ft). At this point, it searches over airspeed 
V (and possibly thrust setting tt ) so that Eq. (A. 39) is minimized. For 
the 727 ’s JT8D engines, thrust is governed by EPR settings which vary 
between 1.1 (idle thrust) and some maximum value less than 2.4. The air- 
speed has an upper limit governed by 

a) . 0.89 Mach structural limits, 

b) . 250 kt (IAS) below 3048 m (10000 ft) for ATC restrictions, 

c) . / 2g(E-h) which insures that the aircraft climbs, and 

d) . the upper value shown in Fig. A, 6 where max thrust 

equals drag. 

The lower limit is governed by 

a) . Vj^, the lower value shown in Fig. A. 6 where max thrust 

equals drag, 

b) . 0.1 Mach 

c) . 5 ft/sec less than the previous value of V to limit large 

jumps in flight path angle. 

The Fibonacci search technique is again used to determine V and tt 
which minimize Eq. (A. 39) for the fixed value of energy E. The value 
chosen for airspeed is accurate to within .0056 Mach, and EPR is accurate 
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to within .009. Associated with these values of V and it are values of 
energy rate E (Eq. (A. 16)) and altitude h; 

h = E -V^/2g . (A. 40) 

From these, approximate values of time, range, flight path angle, and 
fuel burned are obtained from 

At = AE/E, (A. 41) 

sin Y = (Ah/At)/V, 

X = EAx ; Ax = (V cos Y + V^(h)) At 

F = EAf ; Af = w At. 

The above process is repeated by stepping along energy in incre- 
ments of Ae until E^ is reached. The last value of Eq. (A. 39) is 
stored for possible use in evaluating the cruise distance. 

The above climb optimization procedure is repeated with ^ = 1.5^^, 
l.Olii^^, and perhaps other values until the total range of flight converges 
to the appropriate value. This is discussed in further detail later. 


Descent Optimization 


The descent optimization is very similar to the climb optimization 
with regard to the equations which are evaluated. The optimization pro- 
cess requires estimated values of E^ and at the beginning of descent, 
and an estimate of weight at the end of descent. The method used 
to obtain these estimates is discussed in the next section. 


If there is a cruise portion of flight, fuel will be burned off during 
cruise. Thus, the value of E^£> and at the beginning of descent 
will be different than at the beginning of cruise. If there is no cruise 
portion, then these values will be identical. 

The descent profile is obtained by starting at the. final energy 
state and then going backwards in time. The energy rate is constrained 
to be negative with respect to forward time. 
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Similar descent profile constraints exist on airspeed as for those 
of the climb profile. The thrust level is on or near the idle value 
during descent. 

Cruise Fuel Burn 


To estimate the final weight during cruise and landing (W^) , 

the following steps are taken: 


1 ). 

2 ) 


3 ). 


4 ). 


5 ). 


6 ). 


7 ). 


8 ). 

9 ). 


Determine , the initial cruise cost based on the initial 
c 

cruise weight obtained from the climb optimization. 

Use the initial cruise weight and ^ to compute the fuel flow 
rate w(i{j^) . 

Estimate the cruise range d^ by the empirical equations. 


c copt 

d = b.p"'* + b„P^ + b.P^ + b.P + b. • 
cl 2 3 4 5 


Compute the cruise fuel as 


(A. 42) 


F 

c 


= w(ip ) d /(V + V (h )) . 

c c c w c 


(A. 43) 


Estimate the average cruise weight as 

W = W - 0.5F (A. 44) 

c c c 

Use the cruise table to obtain the corresponding cruise cost 
altitude h^, fuel flow rate wOp^), airspeed V^, and wind 
speed V^(h) • 

Recompute Eq. (A. 43) and then find the final cruise weight, 




(A. 45) 


Use the value in the cruise tables to obtain 
As with the climb, set ({^ = 1.5 

Use this value of to obtain h ^ and E from the cruise 

cf cf 

tables. These are the end conditions for the descent trajectory 
obtained backwards in time. 
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10). Estimate the landing weight from the empirical formula 


P = 1.5, (A. 46) 

"f ■ "cf - + '=2’’ + • 


The values of tjj, and obtained by the above procedure are used for 
obtaining the optimum descent trajectory. The descent portion of the 
Hamiltonian is of the form 


C, w + C. 


«dn(E> = 




(V + 




(A. 47) 


this function is also minimized at each of the given values of energy. 


After the first descent profile is completed, a new estimate of 
cruise distance is obtained by using Eq. (A. 35), or 


Then, step (4) above is repeated to obtain an improved cruise fuel burn. 
Then, the Improved landing weight estimate is 


Wf = 


w. 

1 


- (F, 


up 


+ 


F + F , ) 
c dn 


(A. 49) 


The landing trajectory is reoptimized with this new value of 
landing weight. Then, improved values of total range traveled, time 
required, and fuel burned during climb, cruise, and descent are made. 


For short range flight, the above steps assumed that a Type 1 trajec- 
tory is generated because thrust is constrained to maximum value during 
climb and idle value during descent. If thrust is free, then a Type 2 
trajectory will result, with no cruise portion. For this case, the steps 
required to estimate cruise distance d^ and final cruise cost, weight, and 
energy can be eliminated. 
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Cruise Cost Estimation 


The first climb and descent profiles are generated with = 1.5i|^ 
(Ecopt)* next set is generated with = 1.01 i|j • Each of 

these values of cruise cost have an associated range on the curve shown 
in Fig. A. 9. If the total range desired is greater than the value 

obtained using l.Olijj (^^opt^ ’ assumed that the optimum cruise 

altitude and energy are reached. Then a third set of climb and descent 
profiles is generated using In this case, the cruise distance 

is computed so that the desired overall range is exactly achieved. 


If the desired range is between R . and R in Fig. A. 9, then an 
^ min max ^ 

iterative process is used to obtain ^(E^) and the associated desired 
range. Iterations are stopped when the total range traveled is within 
some e of the desired range. (In OPTIM, e is set at 5 n. mi.) 


85 



Flight Range - x, (n.mi.) 


Figure A. 9. Relationship between the Cruise Cost Parameter and 
the Associated Range of Flight 




APPENDIX B 


AIRCRAFT EQUATIONS OF MOTION AND AUTOPILOT MODELS 

The objective of the TRAGEN program is to simulate a 727 aircraft 
being steered to fly along an input reference trajectory. This trajectory 
may either be a climb or descent profile. The simulation must be accurate 
enough such that the performance of the aircraft (in terms of fuel burned 
and time required to reach the destination point) is adequately determined, 
as measured from the output. Adequate accuracy is obtained with a five- 
state variable longitudinal aircraft model. 

The purpose of this appendix is to present the analytical expressions 
upon which the simulation was developed; this is done in two parts. The 
first section below defines the overall system and presents the differential 
equations of motion and fuel burn. The second section describes different 
methods for generating typical guidance commands and autopilot equations. 

Equations of Motion and Fuel Burn 

To examine the vertical profile of the aircraft (i.e., altitude and 
airspeed vs range) , the longitudinal equations of motion are of primary 
importance. The short period equations of motion and the throttle dynamics 
are ignored. Thus, the control variables in this longitudinal plane are 
the angle-of-attack a and the magnitude of the thrust vector T. These 
quantities are shown with respect to aircraft airspeed V , lift L, drag D, 

EL 

weight V, and flight path angle Y in Fig. B.l. 

The kinematic equations of motion of the aircraft in the longitudinal 
plane are 

X = V cos Y + V , (B.l) 

a w 

h = V sin Y > 
a ' ’ 
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Figure B,1 Vector Diagram of Longitudinal Forces 


where 

X - distance, or range, measured on the ground, 

h <- altitude, 

V - wind speed, 

w 


The inertial speed along the airspeed vector V is 

= V + V COSY ' (B.2) 

law' 

From Fig, B.l, the time rate of change of this vector for constant y is 

• * f 1 

= V + V COS Y = — (T cos a - D - W sin y) . (B.3) 

law'm ' 

The time rate of change of the wind speed is 


V 

w 


Substituting Eq. 
V 

a 


3V 

3h ^ ’ 

3V 

w 

V sin "Y 
3h a ' 

(B.4) into Eq. (B.3) and solving for V leaves 

SL 

1 

— (T cos a - D - W sin y) “ V sin y cos y 

m on a 


(B.4) 


(B.5) 
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Also, from Fig. B.l, one can write 

h = — (L cos Y - W + T sin ( y + a ) - D sin Y ), (B,6) 

m 

Equations (B.5) and (B.6) represent the kinetic equations of motion of 
the aircraft. 

The remaining term that must be accounted for is the time-varying 
weight of the aircraft. Specifying the thrust also specifies the fuel 
burn rate w. Thus, the weight changes at the rate 

W = -w . (B.7) 


Equations (B.l), (B,5), (B,6), and (B.7) are the five basic equations 
used to represent the longitudinal dynamics of the aircraft. Lift and 
drag (L and D) are computed as functions of V , h, and a . Fuel flow 

Si 

rate w is a function of V , h, and thrust T. 

a ’ 

Further refinement could be added to these equations to include the 
effects of the following: 

1) . throttle dynamics (Including transient fuel flow rates); 

2) . relationship between throttle position, EPR setting, and thrust; 

3) . short period dynamics relating time rate of change of angle- 

of-attack, pitch rate, and pitch angle to elevator deflection; 

4) . required turning (lateral) motion for flying over fixed waypoints; 

and 

5) . lateral wind and gust effects. 

However, these effects are considered to be of second order, and not re- 
quired for the intent of this simulation. For a more exact autopilot 
simulation, they would be required. 
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The flight path angle is defined as 


Y = sin ^ (h/V ) 

a 

By dif ferentating this expression and using Eqs. (B.5) and (B.6), 

1 2 

Y = -n- (T sin a - W cos Y + L + m V sin y) • 

mV 9h a 

a 

Equation (B.9) can be used in place of Eq. (B,6), 


(B.8) 

one obtains 
(B.9) 


Steering Procedures 


The reference trajectories which are given to be followed consist of 
a sequence of points containing values of time, range, altitude, airspeed, 
flight path angle, specific energy, weight, and other variables. Any of 
these quantities which is measurable and monotonically changing can serve 
as the independent variable. To minimize airborne computer memory re- 
quirements, it is important to make the stored data representing the re- 
ference trajectory as compact as possible. 

In this study, a set of steering equations is used to take points 
from the reference trajectory, convert these points to reference trajectory 
commands, and then use these commands to set values of the control vari- 
ables. This steering process represents a rudimentary form of an autopilot. 

The steering process consists of commanding the thrust T and angle- 
of-attack a values so that the aircraft follows the reference as closely 
as possible. The system that includes this process is depicted by the 
block diagram in Fig. B.2. Note that flying along a reference trajectory 
consists of steering to connect a series of reference points. When a re- 
ference point is reached, new steering commands must be issued so that the 
aircraft will then be guided to the next reference point. 
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1 



Figure B.2. Elements of the Longitudinal Aircraft Model 


To fly along the reference path, an independent variable is first 
chosen. For this study, two different independent variables were chosen - 
altitude for climb and range for descent. Then, the remaining variables - 
primarily airspeed, flight path angle, and altitude (for range as the in- 
dependent variable) - are stored as tabular functions of the chosen in- 
dependent variable. 

Also, it is possible to fly along a nominal path using two approaches: 

1) . An open-loop approach where the thrust vector is directed 

in such a way over the next period that by the end of that 
period the next reference point is reached. 

2) . A closed-loop approach where the aircraft is continually 

steered to a continuously commanded trajectory which connects 
the reference points. 

Both of these approaches were examined for simulation of flying the climb 
profile. However, good results were only obtained with the closed-loop 
approach. 
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A problem with open-loop steering is that it assumes that constant 
or linearly varying controls will cause the end points of a reference 
profile to be connected. This assumption does not account for perturba- 
tions due to wind, etc. along the way. Although the open -loop methods 
produce paths which have roughly correct values of airspeed and altitude 
at given range values, there were large excursions from the .ref erence 
flight path angle for the climb profiles. 

Another problem with the open-loop approaches was that both a and T 
were varied to achieve fixed values of V and h for given range points. 

For optimum climb, thrust is usually set at the maximum value. Thus, 
usually only a remains as a valid control variable. 

Another consideration for implementing the climb profile is that 
there is no reason why a particular cruise condition (altitude, airspeed) 
has to be achieved when a certain range x is reached. Thus, a more logical 
independent variable is altitude, with range allowed to be a free variable. 

For these reasons, a closed-loop steering approach was devised where 

reference values of flight path angle (with respect to the air mass) and 

airspeed are obtained as functions of altitude. (This assumes that 

altitude is monotonically Increasing during climb.) A perturbation control 

law was set up so that variations in a from a reference value a were 

o 

proportional to variations in y and V from their respective command values. 

Si 

Because y and V tend to change linearly with time, they can be 
considered as ramp functions. Thus, the closed-loop controller should be 
considered to be at least a Type 1 system. From Eqs . (B.5) and (B.9), with 

6V - W cos y6y, (B.IO) 

d. 

6V + W sin y6y. 
a 


no wind, the system perturbation equations are 

m<5V = - T sin a 6a - 6a - 

a 6a 6V 

a 

m V 6y = T cos a 6a + 6a + 

Si oCX d V 

a 
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V 



Figure B,3 Control Loops for Perturbation Control of 
Airspeed and Flight Path Angle. 
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where and Y are the commanded values of V and Y* This is the classical 
ac c a 

proportional-plus-integral controller. Gains are chosen to produce the 
desired response for removal of profile errors. 

To generate the continuous commands and Y^ during climb, the computa- 
tions made at each reference point are 


3y _ 

Vl - 

^n 

3h 

h - 

h 

n +1 

n 

3V 

V 

- V 

a 

a n +1 

an 

3h 

h , - 

h 

n +1 

n 


(B.13) 


Then 


- *>n> IC > 


n 


(B.14) 


av 


V = + (h - h ) 37— . 

a a ' n 3h 

c n 


When the flight path angle is very small (during the initial period of 
flight and when the aircraft levels off at 3048m (10000 ft) to gain speed be- 
fore resuming climb), Eqs. (B.14) do not work well. For these cases, it is 
more appropriate to set 


Y^ = 0., 


(B.15) 


V, = 


n +1 


and use the control law 


^4 

6 a = (K 3 + — ) (Y^ - Y) 


(B.16) 
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The above method worked quite well in causing the simulated profile to 
closely follow the reference path. Only one set of gain values was 
sufficient for the entire trajectory. 

For decending flight, the thrust again is usually constrained (idle) 
for optimum performance. Also, for this case, the main concern is to reach 
a fixed altitude when range-to-go to the destination point is a certain 
value. Thus, above 10000 ft, the airspeed can be allowed to be a free 
variable. For this case, only inertial flight path is required to 
be controlled. 

To generate a continuous command the computation made at each 

reference point is 


Ylc " - V ^ <Vl - 

Then, the control law is similar to Eq. (B,16), i.e, 

6 a = (K 3 + — ) (Yj^ - Yj) . (B.18) 

where inertial values of flight path angle are used rather than those with 
respect to the air mass. Equations (B.17) and (B.18) form the basis for 
closed- loop control of descending flight. Again, one set of gains in 
sufficient for the entire descent profile. 
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